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Abstract. The main purpose of the paper is to provide a survey of our recent studies on soliton 
solutions of the Kadomtsev-Petviashvili (KP) equation. The KP equation describes weakly dis- 
persive and small amplitude wave propagation in a quasi-two dimensional framework. Recently a 
large variety of exact soliton solutions of the KP equation has been found and classified. These 
solutions are localized along certain lines in a two-dimensional plane and decay exponentially ev- 
erywhere else, and are called line-solitons. The classification is based on the far-field patterns of 
the solutions which consist of a finite number of line-solitons. Each soliton solution is then defined 
by a point of the totally non-negative Grassmann variety which can be parametrized by a unique 
derangement of the symmetric group of permutations. Our study also includes certain numerical 
stability problems of those soliton solutions. Numerical simulations of the initial value problems 
indicate that certain class of initial waves asymptotically approach to these exact solutions of the 
KP equation. We then discuss an application of our theory to the Mach reflection problem in 
shallow water. This problem describes the resonant interaction of solitary waves appearing in 
' the reflection of an obliquely incident wave onto a vertical wall, and it predicts an extra-ordinary 

four-fold amplification of the wave at the wall. There are several numerical studies confirming the 
• prediction, but all indicate disagreements with the KP theory. Contrary to those previous numer- 

Ch ical studies, we find that the KP theory actually provides an excellent model to describe the Mach 

,1^^ reflection phenomena when the higher order corrections are included to the quasi-two dimensional 

^ approximation. We also present laboratory experiments of the Mach reflection recently carried 

out by Yeh and his colleagues, and show how precisely the KP theory predicts this wave behavior. 
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1. Introduction 



It is a quite well-known story that in August 1834 Sir John Scott Russel observed a large solitary 
wave in a shallow water channel in Scotland. He noted in his first paper (1838) on the subject that 

/ was observing the motion of a boat which was rapidly drawn along a narrow channel 
by a pair of horses, when the boat suddenly stopped - not so the mass of water in 
the channel which it had put in motion; it accumulated round the prow of the vessel 
in a state of violent agitation, then suddenly leaving it behind, rolled forward with 
great velocity, assuming the form of a large solitary elevation, a rounded, smooth and 
well defined heap of water, which continued its course along the channel apparently 
without change of form or diminution of speed 

This solitary wave is now known as an example of a soliton, and is described by a solution of the 
Korteweg-de Vries (KdV) equation. The KdV equation describes one-dimensional wave propagation 
such as beach waves parallel to the coast line or waves in narrow canal, and is obtained in the leading 
order approximation of an asymptotic perturbation theory under the assumptions of weak nonlinear- 
ity (small amplitude) and weak dispersion (long waves). The KdV equation has rich mathematical 
structure including the existence of A^-soliton solutions and the Lax pair for the inverse scattering 
method, and it is a prototype equation of the 1 + 1 dimensional integrable systems. In particular, 
the initial value problem of the KdV equation has been extensively studied by means of the method 
of inverse scattering transform (1ST). It is well known that a general initial data decaying rapidly 
in the spatial variable evolves to a number of individual solitons and weakly dispersive wave trains 
separate from the solitons (see for examples, [Tl [5T1 IMl HI] V 

In 1970, Kadomtsev and Petviashvilli [17J proposed a 2 -I- 1 dimensional dispersive wave equation 
to study the stability of the one-soliton solution of the KdV equation under the influence of weak 
transverse perturbations. This equation is now referred to as the KP equation. It turns out that the 
KP equation has much richer structure than the KdV equation, and might be considered as the most 
fundamental integrable system in the sense that many known integrable systems can be derived as 
special reductions of the so-called KP hierarchy which consists of the KP equation together with 
its infinitely many symmetries. The KP equation can be also represented in the Lax form, that is, 
there exists a pair of linear equations associated with an eigenvalue problem and an evolution of the 
eigenf unction, which enables the method of 1ST. However, unlike the case of the KdV equation, the 
1ST for the KP equation does not seem to provide a practical method of solving the initial value 
problem for initial waves consisting of line-solitons in the far field. 

It is quite important to recognize that the resonant interaction plays a fundamental role in multi- 
dimensional wave phenomenon. The original description of the soliton interaction for the KP equa- 
tion was based on a two-soliton solution found in Hirota bilinear form, which has the shape of 
"X" , describing the intersection of two lines with oblique angle and a phase shift at the intersection 
point. This X-shape solution is referred to as the "0"-type soliton, where "O" stands for original. 
In his study of 1977 on an oblique interaction of two line-solitons. Miles [53] pointed out that the 
0-type solution becomes singular if the angle of the intersection is smaller than certain critical value 
depending on the amplitudes of the solitons. Miles then found that at the critical angle, the two 
line-solitons of the 0-type solution interact resonantly, and a third wave is created to make a "Y- 
shaped" wave form. Indeed, it turns out that such Y-shaped resonant wave forms are exact solutions 
of the KP equation (see also IH]). Miles applied his theory to study the Mach reflection of an inci- 
dent wave onto a vertical wall, and predicted that the third wave, called the Mach stem, created by 
the resonant interaction can reach four-fold amplification of the incidence wave. Several laboratory 
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and numerical experiments attempted to validate his prediction of four-fold amplification, but with 
no definitive success (see for examples UHl US SI] for numerical experiments, and [3S1 [571 SS] for 
laboratory experiments) . 

After the discovery of the resonant phenomena in the KP equation, several numerical and experi- 
mental studies were performed to investigate resonant interactions in other physical two-dimensional 
equations such as the ion-acoustic and shallow water wave equations under the Boussinesq approxi- 
mation (see for examples [T51 Uni HSl ESI HH IH] UHl [33 US])- However, apart from these activities, 
no significant progress has been made in the study of the solution space or real applications of the 
KP equation. It would appear that the general perception was that there were not many new and 
significant results left to be uncovered in the soliton solutions of the KP theory. 

Over the past several years, we have been working on the classification problem of the soliton 
solutions of the KP equation and their applications to shallow water waves. Our studies have 
revealed a large variety of solutions that were totally overlooked in the past [H [521 [SI 13 [HI [HQ ^^"^ 
we found that some of those exact solutions can be applied to study the Mach reflection problem 
[SI [23l Ull [46] . Our numerical study ^20j indicates that the solution to the initial value problem 
of the KP equation with certain class of initial waves associated with the Mach reflection problem 
converges asymptotically to some of these new exact solutions, that is, a separation of dispersive 
radiations from the soliton solution similar to the case of the KdV soliton. 

The main purpose of this paper is to present a survey of our studies on the soliton solutions of 
the KP equation. The paper also presents several results for recent laboratory experiments done by 
Harry Yeh and his colleagues at Oregon State University. 

The paper is organized as follows: 

In Section[2] we present the derivation of the Boussinesq-type equation from the three-dimensional 
Euler equation for the irrotational and incompressible fluid under the assumptions of weak nonlin- 
earity and weak dispersion. The purpose of this section is to give a precise physical meaning to 
those assumptions and to explain the existence of a solitary wave solution in the form of the KdV 
soliton. 

In Section [3| we explain the quasi-two dimensional approximation to derive the KP equation, 
and discuss physical interpretation of the KP soliton in terms of the KdV soliton. In order to 
describe the general soliton solutions, we here introduce the r-function which is expressed by a 
Wronskian determinant for a set of N linearly independent functions {/; : i — 1,...,N}. Each 
function fi is a linear combination of the exponential functions {Ej : j = 1, . . . , M}, where Ej = e^' 
with 9j = kjX + kjy — k^t for some kj G M. The r-function in the Wronskian form was found in 

M 

[551 HOI (see also [11] )• Setting /j = ^ o,ijEj, each solution is parametrized by the N x M 

coefficient matrix A = (aij) of rank N. This representation naturally leads to the notion of the 
Grassmann variety Gr(iV, M), the set of A^-dimensional subspaces given by SpanR{/i : i = 1, . . . , N} 
of R*^ =SpanR{£'j : j ~ and each point of Gr(A^, M) is marked by this yl-matrix 

[3a[21[H]. 

In Section [4| we provide a brief summary of the totally nonnegative (TNN) Grassmann variety, 
denoted by Gr+ {N, M) , which provides the foundation of the classiflcation theorem for the regular 
soliton solutions of the KP equation discussed in the next Section. The main purpose of this section 
is to show that the r-function in the Wronskian determinant form described in Section [3] can be 
identified as a point of the TNN Grassmannian cell. That is, a classification of the regular soliton 
solutions of the KP equation is equivalent to a parametrization of TNN Grassmannian cells [S] [5] [5] . 

In Section [5j we present a classification theorem which states that the r-function identified as a 
point of Gr+ {N, M) generates a soliton solution of the KP equation that has asymptotically M ~ N 
line-solitons for y ^ and N line-solitons for y ^ 0. Moreover, these solutions can be parametrized 



The lower dimensional solutions, called (2, 2)-soliton solutions, have been found by the binary Darboux transfor- 
mation in [5]. 
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by the derangements (the permutations without fixed points) of the symmetric group Sm ■ This type 
of solutions is caUed {M — N, 7V)-sohton solution. The derangements then give a paramctrization of 
the TNN Grassmannian cells, and each derangement is expressed by a unique chord diagram. The 
chord diagram is particularly useful to describe the far-field structure of the corresponding soliton 
solution. (See El El ig.) 

In Section pj we explain all the soliton solutions generated by the r-functions on Gr+(2,4) (see 
also |S]). Some of these solutions are useful to describe the Mach reflection problem in shallow 
water. We show that the A-matrix determine the detailed structure of those solutions, such as the 
asymptotic locations of solitons and local interaction patterns. (See [51 [201 15].) 

In Section [T} we present the numerical study of the KP equation for certain types of initial 
waves. In particular, we consider an initial value problem where the initial wave consists of two 
semi-infinite line-solitons forming a V-shape pattern. Those initial waves were considered in the 
study of the generation of large amplitude waves in shallow water [371 US] • The main result of this 
section is to show that the solutions of this particular initial value problem converge asymptotically 
to some of the exact (2,2)-soliton solutions. These results demonstarte a separation of the (exact) 
soliton solution from dispersive radiation in the manner similar to the KdV case. (See [8l [23l [20].) 

In Section[8l we discuss the Mach reflection problem in terms of the KP solitons which is equivalent 
to Miles' theory (assuming quasi- two dimensionality). We flrst show that the previous numerical 
results (see for examples [13111]), which reported a large discrepancy with the theory, are actually 
in a good agreement with the predictions given by the KP theory. However, here one needs to give 
a proper physical interpretation of the theory when one compares it with the numerical results. 
We also present some laboratory experiments of shallow water waves [JS] [37] . We show that the 
experimental results are all in good agreement with the predictions of the KP theory which can 
describe the evolution of the wave-proflle. Finally we demonstrate that the most complex (2,2)- 
soliton solution associated with the r- function on Gr+(2,4), referred to as T-type solution, can be 
realized in an experiment. 



Sir John Scott Russel continued on to say in his book (1865) that 

This is a most beautiful and extraordinary phenomenon: the first day I saw it was 
the happiest day of my life. Nobody has ever had the good fortune to see it before 
or, at all events, to know what it meant. It is now known as the solitary wave of 
translation. No one before had fancied a solitary wave as a possible thing. 

I hope this survey is successful to convince the readers that the two-dimensional wave pattern gener- 
ated by the soliton solutions of the KP equation is a most beautiful and extraordinary phenomena of 
two-dimensional shallow water waves, and one should have no doubt that observing these patterns 
at a beach will bring the happiest moment of one's life. 



2. Shallow water waves: Basic equations 

Let us start with the physical background of the KP equation: We consider a surface wave on 
water which is assumed to be irrotational and incompressible (see for examples |441 [T]). Then the 
surface wave may be described by the three-dimensional Euler equation. 



A0 = O, 
4>i = 0, 



(2.1) 




for < z < ho + fj, 
at z — 0, 

at z = fj + ho, 
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where </) is the velocity potential with the Laplacian A = 9? + 9~ + 9|, g — 980 cm/sec^ the 
gravitational constant, and /iq the average depth. The first two equations of (2.1) implies, 

4>{x,y,z,i) = cos{z^JX^) ^{x,y,i) with A_L==a?+5~, 
where ip{x,y,t) — (f>{x,y,0,t). In the linear limit, the system can be written in the form, 

cos(/ioyA7) -tpi + 9V = 

This then gives the dispersion relation. 



oj"^ ~ gk tanh/iQ^ = Cg^"'^ I 1 h^k 



(2.2) 



where k yk"^ + k^ and the speed of the surface wave cq — \J gh^ (e.g. cq = 70 cm/sec when 

ho = 5 cm). 

Let us denote the following scales: 

Aq ^ horizontal length scale 

ho vertical length scale = water depth 

oo ~ amplitude scale 

The non-dimensional variables {x, y, z, t, rj, (f)} are defined as 



(2.3) 



X = Aoa;, y = AqJ/, z = Hqz, i = —t, 

Co 



Then the shallow water equation in the non-dimesional form is given by 



(/)^^+/3A_L(/) = 0, 
=0, 

rjt + aV_L(/i • V_L7? = -^02 

where the parameters a and /3 are given by 

ao 



and 



for < z < 1 + ar^, 

at z = 0, 

at z = 1 + ck77. 



The weak nonlinearity implies a ^ 1, and the weak dispersion (or long wave assumption) implies 
/3 <C 1. With a small parameter e ^ 1, we assume 

a-/3 = 0(e). 

As in the previous manner, (j) can be written formally in the form, 

(l){x,y,z,t) = cos ip{x,y,t), 

which leads to the expansion. 
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Then the equations at the surface gives the fohowing system of equations, sometimes called the 
Boussinesq-type equation, 



(2.4) 



/3 



Here we have omitted _l sign. Eliminating in (2.4), we obtain the so-called isotropic Benney-Luke 
equation [5], 



^AV 

3 ^ 



0(e2). 



One can then write this equation in the following form up to the same order. 



(2.5) 



1 - |a ) - A^ + a (VV- • V?At + V • (i/^tV^)) = 0(e2). 



This is a regularized form of the two-dimensional Boussinesq-type equation for the shallow water 
waves. Note here that the dispersion relation of this equation is given by 



L0^ = 



1 



Hl.2 



3 



which agrees with (2.2 1 up to 0{e^). 

It is also well-known that the Boussinesq-type equation can be reduced to the KdV equation for 
a far field with a unidirectional approximation: Let x be the coordinate which is perpendicular to 
the wave crest of a linear shape solitary wave, i.e. 

X = x cos ^I'o + y sin ^fp, 

where (cos ^Po, sin ^"0) is the unit vector in the propagation direction. Then using the far-field 
coordinates in the propgation direction, i.e. 

£.-=X-t, T = et, 



the Boussinesq-type equation (2.5) becomes the KdV equation. 



From the first equation of (2.4), we have 77 = + ^(^)- Then the KdV equation can be expressed 
in the form with physical coordinates. 



(2.6) 



3co cohl ^ 
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The one-soliton solution of the KdV equation is then given by 
(2.7) 



fj = ciQ sech 



'3ao 



X - Co 1 + 



2ho 



0. 



t-Xa 



where oq > and xo are arbitrary constants. One should note that any line-solitary wave in the Euler 
equation (2.1) can be (at least locally) approximated by this soliton under the assumption of weak 
dispersion and weak nonlinearity. This remark will be important when we compare any numerical 
results of the Euler equation or the Boussinesq-type equation with those of the KP equation. 

3. The KP equation 

In this section, we give some basic property of the KP equation, and introduce the soliton solutions 
relevant to shallow water wave problem. In particular, we discuss the physical aspect of the KP 
equation which is derived by a further assumption called quasi-two dimensional approximation (see 
for examples [T71 IT]'). 
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3.1. Quasi- two dimensional approximation and one soliton solution. Let us now assume 
a quasi-two dimensionality with a weak dependence in the y-direction, and we introduce a small 
parameter 7 so that the y-coordinate is scaled as 



(3.1) 



C \/lV, with 7 = 0(e). 



Then the system of equations (2.4) becomes 



Vt + '^xx - a {ipt^x)x - ^4'xxxx + 7V'CC = C'(e^). 
u 

Now we consider a far field expressed with the scaling, 
(3.2) ^ ^ X — t and t — et. 

Then the above equations have the expansions, 



Eliminating 77, we obtain 



Noting T] = + 0{e), we have the KP equation for at the leading order. 



(3.3) 



2e77^ + 3077775 + ^77555 ) + 7?7<;(; =0(6^) 



In terms of physical coordinates, the KP equation is given by 

^77-- 

'Ixx. 



+ y'^TO = CI- 



As a particular solution, we have one line-soliton solution in the form, 

' ' 1 



(3.4) 



f] = ao sech 



' 3ao 
4^ 



X + y tan — co 1 



2/io 



tan"^ W - i(, 



where Og > 0, V&q and xq are arbitrary constants. One should note that the KP equation is derived 
under the assumption of quasi-two dimensionality, that is, the angle ^Pq should be small of order 
0{e), and the solution (3.4) becomes unphysical for the case with a large angle. This can be seen 
explicitly by writing it in the following form in the coordinate perpendicular to the wave crest, i.e. 
X — x cos 5*0 + y sin ^0, 



(3.5) 



77 = ao sech 



3ao 



4/i^ cos2 *n 



X - Co cos *o 1 



2hn 



1 



tan^ ^-0 t-xo 



Noting that cos^fo = 1 — 5 tan^ 4*0 + C(e^) with v^o ~ 0{e), the velocity of the soliton has the 
corrected form up to O(e^), i.e. 

1 



cos ^-o I 1 



2hn 



tan^ *o I = 1 



2ho 



O(e^), 



which does not depend on the angle up to ©(e^). This is consistent with the assumption of the 
quasi-two dimensionality. 

We also note that the line-soliton of (3.4) does not satisfy the KdV equation (2.6 1 except the 
case with = 0. Now comparing the KP soliton (3.4) with the KdV soliton (2.7 1, one can find 
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the correction to the quasi-two dhncnsional approximation, that is, the ampUtudc Og in (3.4) is now 
corrected to 

ao 



(3.6) 



ao 



cos^ * 







With this correction, the KP sohton (3.4) now satisfies the KdV equation (2.6) up to C(e^). This 
correction becomes quite important when we compare our KP resuhs with numerical resuhs of the 
Euler or Boussinesq-type equations, that is, (3.6) gives the relation between the amplitudes of the 
KP soliton and the KdV soliton (see Section [8| for the details). 

3.2. Soliton solutions in the Wronskian determinant. In order to give a general scheme to 
discuss the soliton solutions, we first put the KP equation (3.3) in the standard form, 

(3.7) {Aut + Quux +uxxx)x + ^'^YY ^0, 

where the new variables (X, Y, T) and u are related to the physical ones with 



(3.8) 



2ho 



X — Cot = ho X, 



y = hoY, i^f^T. 

2co 



Hereafter we use the lower case letters (x, y, t) for {X, Y, T) (we do not use the non-dimensional 
variables in (2.3), and the KP variables can be converted to the physical variables directly through 
the relations ( |3.8| )). We write the solution of the KP equation ( |3.7[ ) in the r-function form, 

(3.9) u{x,y,t) = 2dl\nT{x,y,t). 

where the r-function is assumed to be the Wronskian determinant with N functions fi's (see for 
examples gOl Hi [H [H] ) , 



(3.10) 



T = Wr(/l,/2,...,/Ar) 



/l 
/2 



n 



(1) 



N 



f(l) 



JN 



JN 



with /j^"-* = The functions {/j = ,N} satisfy the linear equations dyfi = d^fi and 

dtfi — ^d^fi, and for the soliton solution, we take 



(3.11) 



M 



with Ej = := exp [kjX + k'jy - kjt) . 

of rank N, and we assume the real parameters 



Here A (aij) defines an x M matrix A — (a^j) 
{kj : j — 1, . . . , M} to be ordered, 

ki < k2 < ■ ■ ■ < kM- 

One should emphasize here that we have a parametrization of each soliton solution of the KP equation 
in terms of the /c-parameters and the A-matrix. Then the classification of the soliton solutions is to 
give a complete characterization of the r-function (3.10) with the exponential functions in (3.11). 
This is the main theme in |22l El [H [9] and will be discussed in the following sections, 
r-function, the KP c 

4(rra;t - T^Tt) +TTx 



(3.12) 



In terms of the r-function, the KP equation (3.7) is written in the bilinear form, 

2 



4r.rj, 



3r" 



2,{TTyy - r|^) = 



To show that the r-function (3.10) satisfies this equation, we express the derivatives of the r-function 
using the Young diagram. Let Y be given by the partition Y = (Ai > . . . > A„) where A^'s represent 

n 

the numbers of boxes in Y, and \Y\ denote the total number of boxes, i.e. \Y\ = ^ A^. Denote r in 

1=1 

the A^-tablet, 

r = r0-(O,l,2,...,iV-l), 
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where the numbers describe the orders of the derivative of the column vector (/i, . . . , /n)'^ in the 
T- function (3.10). Then the number of boxes in each row of Y can be found by counting the missing 



numbers which are less than the corresponding number in Ty- For example, T|ip represents 

=(0,l,2,...,7V-3,7V_l,Ar + i). 

For the number + 1, two numbers N — 2 and N are missing, and this gives m . For the number 
N — 1, one number A'^ — 2 is missing and this gives □ . In terms of the determinant, r |ip represents. 



fi 

/2 



{N~3) AN-1) 



N 



J2 



AN-3) 
JN 



/r 



n 



(N+l) 



AN-1) AN+1) 
J2 J2 



AN~1) AN+1) 
JN JN 

_ f(2) 



With this notation and using the equations for fi, i.e. dyfi — fi^',dtfi — — /P\ the derivatives of 



the r-function (3.101 are given by 



Ten = ^{Txx 



Then the bilinear equation ( 3.12[ ) can be written in the form, 
(3.13) '^<I>'^\S +7-mTg =0. 

This equation is nothing but the Laplace expansion of the following 2N x 2N determinant identity. 



3^ 



-Ty) 
rt) 



h 

fN 






n 



{N-2) 



h 



AN~2) „ 
Jn JN 







N 



/i 



{N-S) AN-l) 



{N+l) 



AN-3) AN-1) AN) 

Jn JN JN 



{N-3) AN-1) 



AN+1) 
JN 



/i 



AN-3) 

Jn 



/r 



AN-1) 

Jn 



j.{N) j.{N+l) 



AN) 

Jn 



AN+1) 

Jn 



so that (3.13) is identically satisfied. The relation (3.13) is the simplest example of the Pliicker 



relations (see below), and it can be symbolically written by 

(3.14) C(l, 2)^(3, 4) - ^(1, 3)^(2, 4) + ^(1, 4)^(2, 3) - 0, 

where the Young diagrams are expressed hy Y — (j — 2, i — 1) for the symbol These symbols 

C(jii ■ • ■ iJn) are the so-called Pliicker coordinates of the Granssmann manifold Gr(iV, M). In the 
next section, we outline the basic information for the Grassmannian Gr(iV, M) which will provide 
the foundation of the classification theory for the soliton solutions of the KP equation [22l [8] . 



Example 3.1. Let us express one line-soliton solution (3.4) in our setting. Here we also introduce 



some notations to describe the soliton solutions. The soliton solution (3.4) is obtained by the r- 
function with M = 2 and iV = 1, i.e. 



T = /i = aiiEi + ai2-B2. 



Since the solution u is given by (3.9), one can assume an = 1 and denote ai2 = a > 0. Then 



El + aE2 = 27^65(^1+^=^) cosh i (01 - - Ina) , 
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with the 1x2 A-niatrix of the form A — (1 a). The parameter a in the A-matrix must be a > for 
a non-singular solution and it determines the focation of the sohton solution. Since a — leads to 
a trivial solution, we consider only a > 0. Then the solution u = 29^(lnT) gives 

u = ^(fci - fc2)^sech^ ^(6*1 - 6*2 - Ina). 

Thus the solution is localized along the line 61 — 62 = Ina, hence we call it line-soliton solution. 
We emphasize here that the line-soliton appears at the boundary of two regions where either Ei 
or E2 is the dominant exponential term, and because of this we also call this soliton a [1, 2]-soliton 
solution. In Section 5, we will construct more general line-soliton solutions which separates into a 
number of one-soliton solutions asymptotically as \y\ — > c». We refer to each of these asymptotic 
line-solitons as the [i, j]-soliton. The [i, j]-soliton solution with i < j has the same (local) structure 
as the one-soliton solution, and can be described as follows 



sech^ ^ (^['j] ■ ^ ^ + '^[i, 



with some constant ©j'^ jj- The amplitude ^[ij], the wave- vector K[j j] and the frequency ^[ij] are 
defined by 

The direction of the wave- vector K[j j] = (if , ii'^ ) is measured in the counterclockwise sense 
from the j/-axis, and it is given by 



= tan^Pij j] = ki + kj, 



that is, ^'[i.j] gives the angle between the line K^j jj • x — const and the y-axis. Then one line-soliton 
can be written in the form with three parameters ^'[ij] £^nd Xj*. ^j. 



(3.15) u = sech^ y (^x + y tan - q, jjt - xf,^^-] 

with q.j j] = fc| -I- fci/cj + = h^lij] + l^^'^^^lij]- Figure |3.ll we illustrate one line-soliton 
solution of [«, j]-type. In the right panel of this figure, we show a chord diagram which represents 
this soliton solution. Here the chord diagram indicates the permutation of the dominant exponential 
terms Ei and Ej in the r- function, that is, with the ordering ki < kj, Ei dominates in a:; ^ 0, while 
Ej dominates in cc ^ (see Sect ion [4| for the precise definition of the chord diagram) . 

For each soliton solution of (3.15), the wave vector K[jj] and the frequency ^[ij] satisfy the 
soliton-dispersion relation (see (2^)), 

(3-16) m[,,^Kl^^^iKl^^)^ + 3iKl^^r. 

The soliton velocity V[j jj is along the direction of the wave-vector K[,j j], and is defined by ■ 
^[i.i] = ^[i.i]' which yields 

r2r,; A] k- -|- kikj -\- k- 

^^-^'=pcj^'--'^ i+(..:fc,)^ 

Note in particular that since C[i.j] ~ kf + kikj + A;| > 0, the x-component of the soliton velocity is 
always positive, i.e., any soliton propagates in the positive x-direction. In the physical coordinates 
(see (3.8)), this implies that soliton propagates in super-sonic (i.e. the speed of soliton is faster than 
Co = \/gho, because of its nonlinear effect with fj > 0, see Section [2]). On the other hand, one should 
note that any small perturbation propagates in the negative x-direction, i.e., the x-component of the 
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Figure 3.1. One line-soliton solution of [«, j]-type and the corresponding chord 
diagram. The ampUtude ^[ij] and the angle '^[ij] are given by ^[ij] = ^{ki — kj)'^ 
and tan "^[ij] = ki + kj. The upper oriented chord represents the part of [i, j]-soliton 
for y ^ and the lower one for y <C 0. 



group velocity is always negative. This can be seen from the dispersion relation of the linearized KP 



equation for a plane wave 



exp(zk • X — iuit) with the wave- vector k = [k^, ky) and the frequency 



- ^kl 



kl 



from which the group velocity of the wave is given by 



V = Vcj = 



dkx 



du! 
dky 




3 ky 



Physically, this means that the radiations disperse with sub-sonic speeds. This is similar to the case 
of the KdV equation, and we expect that asymptotically, the soliton separates from small radiations. 
We further discuss this issue in Section [7] where we numerically observe the separation. 



Remark 3.2. In the formulas (3.11 1, if we include the higher times t„ in the exponential functions, 
i.e. 



exp 



kj tr 



then the r-function (3.10) gives a solution of the KP hierarchy. The equation for the i„-flow is a 



symmetry of the KP equation, and the r-function with those higher times also satisfies the other 
Pliicker relations which are expressed with the Young diagrams having larger numbers of boxes [30) . 

4. Totally nonnegative Grassmannian Gr+(7V, M) 
In the previous section, we considered a class of solutions which are expressed by the r-functions 



( 3.10 ) with the exponential functions ( 3.11 ). Those solutions are determined by the fc-parameters and 
the ^-matrix. Fixing the /c-parameters, we have a set of M exponentials {Ej — e^^ : j — 1, . . . , M} 
which spans M^^. Then the set of functions {fi : i = 1, . . . ,N} of (3.11) defines an TV-dimensional 



subspace of M*'. This leads naturally to the notion of Grassmannian Gr(iV, M), the set of all A^- 
dimensional subspaces in M.^^ , and each point of Gr(A'^, M) can be parametrized by the A- matrix 
in (3.11). Here we give a brief review of the Grassmann manifold Gr(A^, M), in particular, we 



describe the totally non-negative part of Gr(iV, M). The main purpose of this section is to explain 
a mathematical background of regular soliton solutions of the KP equation. 

4.1. The Grassmannian Gr(A^, M). Recall that the set of the functions fi spans an iV-dimensional 
subspace which is parametrized by an N x M matrix A of rank N , i.e. 



(/i, /2, ■ • ■ , In) — iEi,E2, 



.E 
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Since other set of functions (51, • • • , A'at) = {fi, . . . , fN)H for some H eGLjv(M) gives the same 
subspace, the A-matrix can be canonically chosen in the reduced row echelon form (RREF). This 
then gives an explicit definition of the Grassmannian, 

Gr(7V,M) = GLAr(M)\7W 

where AiNxM(^) denotes the set of N x M matrices of rank N. The canonical form of A is 
distinguished by a set of pivot columns labeled by / = {ii, Z2, . . . , zat}, 1 < ii < i2 < ■ ■ ■ < iN ^ M 
such that the N x N sub-matrix Aj formed by the column set I is the identity matrix. Each N x M 
matrix A in RREF uniquely determines an A/'-dimensional subspace, thus providing a coordinate 
for a point of Gr(A'', M). The set Wj of all points in Gr(A^, M) represented by RREF matrices 
A which have the same pivot set / is called a Schubert cell which gives the decomposition of the 
Grassmannian, the Schubert decomposition, 

(4.1) Gt{N,M)= y Wi, I^{ii,i2,...,iN}. 

l<ii<i2<-..iN<M 

For example, if I = {1,2,..., N}, then the Schubert cell Wj contains all A matrices whose RREF 
is given by 

/ 1 ••• * ••• * \ 
1 ••• * ••• * 

\0 ••• 1 * ••• */ 
where the N(M — N) entries of the right-hand block are arbitrary real numbers. This particular 
Schubert cell is often referred to as the top cell which has the maximum number of free parameters 
marked by *. It follows from this that the dimension of Gr(7V, M) is N{M — N). The number of 
free parameters for an A-matrix in RREF with a given pivot set / = {ii, . . . , in} is equal to the the 
dimension of the cell Wj, and is given by 

TV 

dim = N{M - N) ~ T,{in-n). 

n=l 

Note here that the index set I = {ii, . . . ,ijv} can be expressed by the Young diagram with Yj = 
{iN — iV, . . . , «2 — 2, «! — 1), and then codim Wi = \Y\. 

Example 4.1. The Schubert decomposition of Gr(2,4) has the form, 

Gr(2,4)= □ Wiij}. 

l<i<j<4 

There are six cells Wj with dim Wj = 7 — {i+ j), and are listed below: 

(-)^n.>=^«={(; ; : :)}' i^)w,.,,=w.=[(i ; ; ;)}, 
(c);^,i,>=^-={(; ; I ;)}> i<^)w,.,sy=w,={(^^ ; ; :)}, 

The top cell M^{i,2} has four free parameters which gives dim Gr(2, 4), while the bottom cell ^{3^4} 
is 0-dimensional and corresponds to a single point of the Grassmannian. 

We also note that each cell in the Schubert decomposition can be parametrized by a unique 
element of Sm, the symmetric group of permutations for M letters. The group Sm is generated by 
the adjacent transpositions Sj := + 1), i.e. 

Sm = {si,S2, • . . , sm-i), 
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with sf = e, the identity element, SiSj = SjSi if \i — j| > 1 and (5^5,^+1)'^ = e. Let Pjv be a maximal 
parabolic subgroup of Sm generated by sj's without the element sm-n, i-e. 

Pn ■= (si, . . . , sj\/-Af-i, sm-n+1, ■ ■ ■ , sm-i) — Sm-n x Sn- 

Then the pivot set / = {ii, 22, ... , i^v} parametrizing the Schubert cell Wj can be uniquely labeled 
by a minimal length representative of the coset , 

S'l^'' = {the reduced words ending with sm-n }■ 

Pn 

Namely, we have the Schubert decomposition of Gr(iV, M) in terms of the coset S^j^\ 



Gr(7V,M)= y W^. 



(N) 



ires;, 

where the dimension of the cell is given by the length of the permutation, i.e. dim W^r = £{7t) 
For example, in the case of Gr(l,3), we have Sg^"* = {si, S2) / {si) = {e, S2, S1S2}, 

1 2 0\ J. A [2] 3\ ^ /[T] 2 3 
123/ ^132y 1^312 

Here |T] represents a pivot, so that we have 

T^e^ {(0,0,1)}, ly,, ={(0,1,*)}, 1^,,,, ={(1, *,*)}. 

Also in the case of Gr(2,3), we have ^3 ' — {si, S2) / {S2) = {e, si, S2S1}, 

1 [2}{3]\ ^ (^0 2 0\ ^ /^E 3 
123/ ^213/ V23I 

and the Schubert cells W^^ are given by 

Note in particular that the last elements in the above examples have no fixed points, and they are 
called derangements. As we will show that each derangement of Sm parametrizes a unique line- 
soliton solution generated by the r- function of the form (3.10). It is important for our purposes to 
remark that each permutation tt S Sm with marked pivot positions can be uniquely expressed by the 
chord diagram. This permutation is the decorated permutation defined in |38j for a parametrization 
of the totally non-negative Grassmann cells. 

Definition 4.2. A chord diagram associated with tt e Sm is defined as follows: Consider a line 
segment with M marked points by the numbers {1,2,..., Af } in the increasing order from the left. 

(a) If i < 7r(i) (excedance), then draw a chord joining i and 7r(i) on the upper part of the line. 

(b) If J < 7r(j) (deficiency), then draw a chord joining j and 7r(j) on the lower part of the line. 

(c) If I = Tr{l) (fixed point), then 

(i) if Z is a pivot, then draw a loop on the upper part of the line at this point. 

(ii) if Z is a non-pivot, then draw a loop on the lower part of the line at this point. 

The dimension of each Schubert cell of Gr(A^, M) can be also found from the chord diagram, and it 
is given by 

dimV^V —N + {# of crossings} + {# of cusps in the lower part} 
— {# of loops in the upper part}. 

Here we say that the point marked by j is a "cusp" , if 7r(j') < j — n{k) < fc or fc < 7r(/c) — j < for 
some k. In particular, the point j is a cusp in the lower part of the diagram, if < j = n{k) < k 
(see [inillS]). 
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Example 4.3. Consider the case of Gr(2,4). The Schubert cells are marked by the pivots 

with 1 < * < J < 4, and the permutation representations are given by 

1 2 a A @ 3 0\ ^ (^0 2 3 

123 4/ ^1324/ ^3124 

S3 i S3 i 

1 @ S 4\ ^ /0 2 4' 
13 42/ ^314 2 

S2 i 

@ 3 4' 
3 4 12 

The chord diagrams are shown below, and the points with filled circle indicate the pivots for those 
cells: One should note here that each fixed point corresponds to a loop of the diagrams, and the 



1 k2 { 


) 


1 ( 


j k3 k4 







diagrams without loops are associated with the derangements of the permutation group. 

As we will show, each chord (not loop) identifies a line-soliton for y ^ (or y <C 0) corresponding 
to the location of the chord in the upper (or lower) part of the chord diagram. For example, in the 
case TT = (3 ^ 4 2)' '^^ have [1,3]- and [3, 4]-soHtons in y and [1,2]- and [2, 4]-solitons in y <C 0. 

4.2. The Pliicker coordinates and total non-negativity. We here describe the totally non- 
negative (TNN) Grassmannian Or"*" (iV,M) as a subspace of Gr(iV, M). Then we will show that 
the T- function associated with Gr"*" (iV, M) is necessary and sufficient conditions for the solution 
generated by the r-function to be regular. 

We first note that the coordinates of Gr(iV, M) is given by the Pliicker embedding into the 
projectivization of the wedge product space A^M*^, i.e. 

Gr(iV,M) ^ P(A^M*^), 

which maps each frame given by [/i, . . . , /jv] GGr(iV, M) to the point on P(A^M*'^), i.e. 

(4.2) /iA-.-A/at^ E aJi,---,JN)E,,A---AEj,. 

l<jl<:.<jN<M 

Here the coefficients £,{ji, ■ ■ ■ , jw) are the N x N minors of the ^-matrix defined by 

CUi,- ■ -Jn) ■■= det|(a„j„)i<„<Ar|. 

Note here that ^(ji , . . . , jV) = 1 where the set {ji , . . . , jat} is the pivot set. Those minors xi{ji , . . . , j^r) 
are called the Pliicker coordinates, which give a coordinate system for the linear space A^M*^ with 
the basis, 

{E,, A • • • A : 1 < ji < • • • < JAT < A/}. 
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Then the Grassmannian structure is determined by certain relations on the Pliicker coordinates, 
called the Pliicker relations: for any two index sets {ai, . . . ,ajv-i} and . . . ,/?jv+i} with 1 < 
ai, Pj < M, they are given by 



(4.3) 



N+l 

E e(«i, 



,aN-i,l3j)£,{l3i, ■ ■ ■ ,l3j, ... ,I3n+i) = 0, 



0. 



where $j implies the deletion of Pj. The Pliicker relations can be derived using elementary linear 
algebra from the Laplace expansion of the following 2N x 2N determinant formed by the columns 
Ai of the matrix A, i.e. 

Aai ■ ■ ■ ^Qjv-i ^/3i • • • ^^jv+l 

••• A^, ■■■ A/3„+, 
The Pliicker coordinates, modulo the Pliicker relations, give the correct dimension of Gr(iV, M) 
which is typically less than the dimension of P(/\ M*^). 

Example 4.4. For Gr(2,4), the Pliicker coordinates are given by the maximal minors, 
^(1,2), ^(1,3), e(l,4), e(2,3), ^(2,4), ^(3,4). 



Taking {ai} = {1}, /32, /^s} = {2,3,4} in (4.3) gives the only Pliicker relation in this case, 

2)^(3, 4) - 3)C(2, 4) + 4)^(2, 3) = 0, 

which is the same as (3.14). Since dim(/\^ M^) — 6, and the projectivization gives dim(P(/\^ M^)) = 
6 — 1 = 5. Then with one Pliicker relation, the dimension of Gr(2,4) turns out be 4, which is 
consistent with the dimension of the top cell W{i 2} as shown in Example 3.2 (case (a)). 



Since each point of Gr(A^, M) is expressed by (4.2), the TNN Grassmannian Gr+(A^, M) is defined 
by the set of iV x M matrices of rank N whose minors, the Pliicker coordinates, are all non-negative, 
i.e. 

Gr+(iV, M) := {A € Gr(7V, M) : ^(ji, . . . ,iiv) > 0, V 1 < ji < • • • < jm < M} . 
Then the most interesting question is to find a parametrization of all the cells in Gr+ {N, M) . 
This question has been solved by Postnikov and his colleagues (see [351111]), and our classification 
theorem of the soliton solutions provides an alternative proof based on a simple asymptotic analysis 
as described in Section [s] (see also [SI [5]). 

4.3. The r-function as a point on Gr+(iV, Af). Expanding the r-function in the Wronskian 
determinant (3.10) by Binet-Cauchy formula, we have 

(4.4) 



Wr(/i,/2,...,/jv) 



E C(«l,*2,. 
l<ii<---<ijv<Af 



■ ,iN)E{ii,i2, ■ ■ ■ ,iN) : 



where ^(11,^2, 
and E{ii, 12, . 
A 



. ,iN) are the Pliicker coordinates given by the maximal minors of the A-matrix, 

, ipf) can be identified with 



E,, A E, 



■■Jn) = Wr(£'jj,£;j2, . . . ,£'j„). Here each E{ii,i2, 
■ Eij^ which is the basis element for i.e. 



N 



Spanj5 {E{ii,i2, 



,tN) 



1 < ii < i2 < ■ ■ ■ < iN < M} = f\\ 



Note here that the sum fcj^ + k^^ 



.. + ki^ should be distinct for distinct sets {zi, 12, . . . , iat} 
in order for the functions {i?(«i,i2, . . . , Jat)} to be linearly independent. It is then clear that the 
T- function given by the Wronskian determinant (3.10) can be identified with a point of Gr(A'^, M), 
and the Wronskian map Wr: [/i,...,/^?] ^ t gives the Pliicker embedding. With the ordering 
ki < ■■■ < k^f, the Wronskian Wr{Ei^, . . . , Ei^) > for ^ < • • • < ijv- Then r eGr+{N,M) 
implies that r-function is positive definite and the solution u{x,y,t) = 29^(lnT) is regular for all 
{x,y,t) G R'^. In order to prove a converse of this statement, we first show the following: Let 
(ti, t2, . . . , tjv/) be the higer times for the KP equation (see Remark |3.2[ and here the first three 
times ti = x,t2 = y and = ~t give the KP variables). 
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Proposition 4.1. Suppose that the r-function is regular for all {ti,t2, ■■■ ,tM) ■ Thenr GGr+(7V,Af). 
Proof. Let us first write the exponential terms, 

/ M \ - 

= exp i J2 kpn + O^j =■■ Eje"^ for j = 1, . . . , M, 

with 9'j e M, i.e. the shifts of i„'s in the exponential functions. Because the fc-parameters are all 
distinct, one can take the coordinates {9i, . . . ,0^j) instead of {ti, . . . ,tM)- The r-function is then 
given by 

N ^ 

l<Ji<---<J]v<M 1^^-^ 

where E{ji , . . . , jjv) = WT{Ej- , . . . , Ej^^ ) . Then one can choose the parameters {9j^ , . . . , 9'j^ ) so that 

N N 

E0l » Eol 

k=l k=l 

for any other choice of the parameters (0;° , . . . ,0'^^). This means that the exponential term having 
this index set {ji, . . . , Jat} is the dominant one in the r-function, while all other parameters are 
of 0(1). Suppose the minor . . . , jat) associated with this index set is negative, that is, r ss 
. . . , jpf)E{ji, . . . , Jn) < 0. Now note that for a; <C 0, the dominant exponential in the r- 
function is E(ei, . . . , cn) > with the pivot set {ei, . . . , ejv} and the ordering A;i < • • • < /c^/, so 
that T « E(ei, . . . , e^v) > 0. This implies that the r-function vanishes at some point in (a;, 2/)-plane, 
and therefore the solution u{x, y, t) is not regular. □ 

It is then clear from the proof that the total non-negativity is not only sufficient but necessary 
for the regularity of the solution. Namely we have the following. 



Corollary 4.1. The solution of the KP equation generated by the t -function in the form (3.10) with 



(3.11) is non- singular for any initial data if and only if t €Gr'^{N,M). 



Thus the classification of the regular soliton solutions is equivalent to a study of the totally 
non-negative Grassmannian. 

Remark 4.5. Since each r-function can be identified as a point on Gr(A^, Af), one can define a 
moment map, fi : Gr(iV, M) — > f)^ 

E ie(ji, . . . ,jN)Ein, . . . ,jN)\HLn + • • • + 

E \aji,---,jN)E{j^,...,jNW 

I<j<---<JN<M 

where Lj are the weights of the standard representation of SL(A/), and f)g is the real part of the 
dual of the Cartan subalgebra of sl(Af ) defined by 

(,^ = SpanR|Li,...,LM : E = o| ^ M^^^^. 

Then the closure of the image of the moment map is a convex polytope whose vertices are the fixed 
points of the Sm orbit, that is, the dominant exponentials. 

5. Classification of soliton solutions 



In this section, we now show the asymptotic behavior of the r-function in (4.4) and then present a 
classification scheme for the regular line-soliton solutions of KP based on the r-function asymptotics 
(see also [S]). 
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/ r|i 




[i,3r><f 




r|2 






\[2,3] 




/ [1,2] 





c 



Figure 5.1. The graphs of 77j(c) = — c). Each [i, j] represents the exchange 

of the order between r\i and r]j. The parameters are given by (fci,/c2,fc3) = 

f_5 _1 3n 



5.1. Asymptotic line-solitons. The r-function of (4.4 1 is given exphcitly by the sum of exponen- 
tial terms with the Wronskians, 

(5.1) T{x,y,t)^ J2 ^{mi,...,mN)E{mi,...,m]y), 

1 < m 1 < ■ ■ ■ < m iv ^ 

where E{mi, . . .,171^) = Wr{E^^,. . with = e"", and 6'™(a;,?/,t) = fc^a; + fc^y - k^t. 

Since u — 29^(lnT), the regularity condition on the hne-sohton solutions requires that the r-function 
does not vanish for all values of x, y and t. To ensure that it is sign-definite, the following necessary 



and sufficient conditions are imposed on the r-function in (5.1). 

(i) The parameters fci, /c2, • ■ • are ordered as fci < ^2 < . . . < km, and the sums ki + kj are 
all distinct. 

(ii) A is a TNN matrix, that is, all its maximal minors are ^(mi, . . . , niN) > 0. 

The asymptotic spatial structure of the solution u{x, y, t) is determined from the consideration of 
dominant exponentials E{mi, . . . ,m7v) in the r-function at different regions of the (x, y)-plane for 
large \y\. The solution u — 29^(lnr) is localized at the boundaries of two distinct regions where a 



balance exists between two dominant exponentials in the r-function (5.11, whereas the solution is 
exponentially small in the interior of each of these regions where only one exponential E{mi, . . . , tojv) 
with a specific index set {mi, . . . , tun} is dominant. Before discussing a general theorem, let us first 
consider the following simple examples which illustrate the resonant interactions among the line- 
solitons. As discussed in Introduction, the resonant interaction is one of the most important features 
of the KP equation (see e.g. [251 [5^ [T5] ). 

Example 5.1. We consider the case with iV = 1 and M — 3, where the 1x3 coefficient matrix A 
is given by 

^ = (1 a b). 

The parameters a, b in the matrix are positive constants; meaning that this A-matrix marks a point 
on Gr+(1, 3), and the positivity implies the regularity of the KP solution. The r-function is simply 
given by 

T^f^Ei+aE2 + bE3. 

Now let us determine the dominant exponentials and analyze the structure of the solution in the 
xy-plane. First we consider the function / along the line x = —cy with c = tan ^' where ^ is the 



angle measured counterclockwise from the j/-axis (see Figure 3.1 1. Then along x = —cy, we have the 
exponential function Ej — exjp[rij{c)y — k^t] with 

(5.2) ^i(c) = kj{kj — c). 



It is then seen from Figure 



5.1 



that for ?/ S> and a fixed t, the exponential term Ei dominates when 
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c is large positive (^P « |, i.e. x — > — oo). Decreasing the value of c (rotating the line clockwise), 
the dominant term changes to E^. Thus we have 

„ , . f fci as a; — >■ — cx3, 
w := Ox m / — > < , 

[ fcs as X — > oo. 

The transition of the dominant exponentials Ei — > E^ is characterized by the condition rji = 773, 
which corresponds the direction parameter value c = tan^E'j]^ 3] — ki + k^. In the neighborhood of 
this line, the function / can be approximated as 

f^Ei+ bEs, 

which implies that there exists a [1, 3]-soliton for y ^ 0. The constant b can be used to choose a 
specific location of this soliton. 

Next consider the case of y 0. The dominant exponential corresponds to the least value of 
Tjj for any given value of c. For large positive c (5" « ^, i.e. x — >■ 00), E3 is the dominant term. 
Decreasing the value of c (rotating the line x — —cy clockwise), the dominant term changes to E2 
when k2 + k^ > c > ki + k2, and Ei becomes dominant for c < ki + k2. Hence, we have for y 

{/si as a: — > —00, 
k2 for -{ki + k2)y <x < -{k2 + ks)y, 
k^ as X — 00. 

In the neighborhood of the line x + (fci + k2)y ^constant, 

which corresponds to a [1, 2]-soliton and its location is fixed by the constant a. The solution in 
?/ <C also consists of a [2, 3]-soliton in the neighborhood of the line x + (fc2 + kz)y =constant, and 
whose location is determined by the locations of other line-solitons. Therefore, we need only two 
parameters a, b (besides the fc-parameters) to specify the solution uniquely, and those parameters 
fix the locations of line x + y tan j] — C^i jT^t = 2:° for [i, j]-soliton. For [z, j] — [1,3] and [2, 3] in 
X ^ 0, we have 

The shape of solution generated by / = i^i + aE2 + bEj, with a = = 1 (i.e. &l t ~ Q three 



line-solitons meet at the origin) is illustrated via the contour plot in Figure 5.2 In this Figure 



one can see that the line-soliton in y ^ labeled by [1,3], is localized along the line 61 = 6^ with 
direction parameter c = fci + fc3; two other line-solitons in y ^ labeled by [1,2] and [2,3] are 
localized respectively, along the phase transition lines with c — ki + k2 and c = k2 + k^. This 
solution represents a resonant solution of three line-solitons. The resonant condition among those 
three line-solitons is given by 

K[i^3] = K[i^2] + K[2,3], ri[i^3] = il[i^2] + ^^[2,3], 

which are trivially satisfied with K[i = {kj — fci,/c| — kf) and = — kf. The resonant 

condition may be symbolically written as 

[1,3] = [1,2] + [2,3]. 

One can also represent this line-soliton solution by a permutation of three indices: {1,2,3} which 
is illustrated by a (linear) chord diagram shown below. Here, the upper chord represents the [1,3]- 
soliton in y ^ and the lower two chords represent [1, 2] and [2, 3]-solitons in y ^ 0. Following the 
arrows in the chord diagram, one recovers the permutation, 

'1 2 3^ 



3 1 2 



or simply n — (312). 



In general, each line-soliton solution of the KP equation can be parametrized by a unique permutation 
corresponding to a chord diagram (see the next subsection). 
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Figure 5.2. Example of (2, l)-solitoii solution and the chord diagram. The k- 
parameters are chosen as (fci,fc2,fc3) = (~ti~ii|)- The right panel is the cor- 
responding chord diagrams. We take the A-matrix A = (1 1 1) so that at i = 
three line-solitons meet at the origin. Each E{j) with j = 1, 2 or 3 indicates the 
dominant exponential term Ej in that region. The boundaries of any two adjacent 
regions give the line-solitons indicating the transition of the dominant terms Ej. 

and the line-solitons are 
Here a = 6 = 1 



The fc-parameters are the same as those in Figure 5.1 



5.1 



determined from the intersection points of the rjj (c) 's in Figure 
(i.e. T = El + E2 + E-s) so that the three solitons meet at the origin at i = 0. 



The results described in this example can be easily extended to the general case where / has 
arbitrary number of exponential terms (see also [26l l4j). 

Proposition 5.1. If f — aiEi + a2i?2 + • • • + umEm with aj > for j — 1,2,..., M, then the 
solution u consists of M ~ 1 line-solitons for ?/ <C and one line-soliton for y ^ 0. 

Such solutions are referred to as the {M — 1, l)-soliton solutions; meaning that {M — 1) line- 
solitons for y ^ and one line-soliton for y ^ 0. Note that the line-soliton for j/ ^ is labeled 
by [1, M], whereas the other line-solitons in ?/ ^ are labeled by [k, k + 1] for k = 1, 2, . . . , M — 1, 
counterclockwise from the negative to the positive x-axis, i.e. increasing Vl/ from — | to As in 
the previous examples one can set ai = 1 without any loss of generality, then the remaining M — 1 
parameters 02, ■ ■ ■ , om determine the locations of the M line-solitons. Also note that the xy-plane 
is divided into M sectors for the asymptotic region with x"^ + y"^ 0, and the boundaries of those 
sectors are given by the asymptotic line-solitons. This feature is common even for the general case. 



Figure 5.3 illustrates the case for a (3, l)-soliton solution with f = Ei + E2 + E^ + E4. The chord 



diagram for this solution represents the permutation tt = (4123) = S4S3S2 € 

Example 5.2. Let us now consider the case with N = 2 and M = 3: We take the A-matrix in 



(3.11 ) of the form, 

1 -b 



A 



Ola 



where a and b are positive constants, that is, A marks a point on Gr+(2, 3). Then the r-function is 
given by 

r = E{1, 2) -I- aE{l, 3) + bE{2, 3). 
In order to carry out the asymptotic analysis in this case one needs to consider the sum of two rjj (c) , 



i.e. r]i,j — rji + rjj for 1 < i < j < 3. This can still be done using Figure 5.1 but a more effective 



way is described below (see the graph of r]{k,c) in Figure 5.5) 



For y ^ 0, the transitions of the dominant exponentials are given by following scheme: 

E{1,2) ^ E{1,3) E{2,3), 
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Figure 5.3. The time evolution of a (3, l)-soliton solution and the corresponding 
chord diagram. The upper chord represents the [1, 4]-soliton, and the lower ones 
represent [1,2]-, [2,3]- and [3, 4]-solitons. The chord diagram shows tt = (4123). 
The intermediate solitons are [1, 3]-soliton a.t t ^ —8 and [2,4]-soliton a.t t — 8, 
respectively. These solitons appear as the resonant (2, l)-type solutions. 



as c varies from large positive (i.e. x — )■ —oo) to large negative values (i.e. x — >■ oo). The boundary 
between the regions with the dominant exponentials E{1,2) and i?(l,3) defines the [2, 3]-soliton 
solution since here the r-function can be approximated as 

r « Eil, 2) + aEil, 3) = 2(^2 - fci)e'^ + 5(»^+«3-e23) cosh i(02 - ^3 + ^23), 
so that we have 

u = 292 Inr w i(/c2 - fe3)'sech"i(02 - ^3 + ^23) , 
where ^23 is related to the parameter of the A-matrix (see below). A similar computation as above 
near the transition boundary of the dominant exponentials E{1,3) and £'(2,3) yields 

T w 2(fc3 - fci)ae*='+5(«i+»2-»i2) cosh i(6ii -62 + 6*12). 
The phases 6*12 and ^23 are related to the parameters of the ^-matrix. 



ks - ki K3 - K2 

For y <ti 0, there is only one transition, namely 

£(2, 3) ^£(1,2), 

as c varies from large positive value (i.e. x — >■ 00) to large negative value (i.e. x — > —00). In this 
case, a [1, 3]-soliton is formed for y ^ at the boundary of the dominant exponentials £(2,3) and 
E{1, 2). The contour plot of the line-soliton solution is shown in Figure 5.4 Notice that this figure 
can be obtained from Figure 5.2 by changing {x,y) — ^ (— x, —y). This solution can be represented 
by the chord diagram corresponding to the permutation tt — (231) shown below. Note that this 
diagram is the n-rotation of the chord diagram in Example |5.1| whose permutation tt = (312) is the 
inverse of tt = (231). 

As shown in those examples, it is now clear that each line-soliton appears as a boundary of two 
dominant exponentials, and with the condition that ki + kj are all distinct for i ^ j, we have the 
following Proposition: 
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Figure 5.4. Example of (1, 2)-soliton solution and the chord diagram. The k- 
parameters are the same as the (2, l)-soliton in the previous figure. The right 
panels are the corresponding chord diagrams. The parameters in the A-matrices 
are chosen as a = ^ and 6 = 1 so that at t = three line-solitons meet at the origin. 



Proposition 5.2. Two dominant exponentials of the t -function in adjacent regions of the xy-plane 
are of the form E{i, m2, . . . , tojv) and E{j, m2, ■ ■ ■ , ra^) for some N—1 common indices TO2, . . . , mpf. 



As a consequence of Proposition |5.2[ the KP solution behaves asymptotically like a single line- 
soliton 

(5.3) u{x, y,t)~\ (fc, - hfsech" \ {9, - 9, + %) , 

in the neighborhood of the line x + (ki + kj)y — constant, which forms the boundary between the 



regions of dominant exponentials i?(z, m2, . . . , totv) and E{j,m2, ■ ■ ■ -,^1^)- Equation (5.3) defines 
an asymptotic line-soliton, i.e. [z, j]-soliton, as a result of those two dominant exponentials. In 
order to identify the set of asymptotic line-solitons associated with a given solution, we need to 
determine which exponential terms E(mi, m2, . . . , mjv) are actually dominant along each line : 
X = —{ki -f kj)y as \y\ — > oo. For this purpose, first note that along a line x = —cy each exponential 
term E{mi, m2, . . . , m^) has the form, 

E{mi,m2, ■ . -jniN) oc exp ^m„(c)y^ , ?7m(c) = k„i{km - c) . 

Thus for y (or <C 0) , the dominant exponential corresponds to the largest (or least) value of the 
sum of 77m„ (c) for each c. When two dominant exponentials E{i, m2, ■ ■ ■ , rnpf) and E{j, m2, . . . , tun) 
are in balance along the direction of the [i, j]-soliton, we have r]i(c) — rij{c) which implies that 
c = ki + kj. Since rjmic) — rii{c) ~ {km — ki){km + h — c) and the fc-parameters are ordered as 
ki < k2 < ■ ■ ■ < kM, we have the following order relations among the other rym(c)'s along c = ki + kj, 



(5.4) 



i]i = Vj < Vrn if m < i or j < m, 
rji = T]j > rjm if i <m < j. 

The relations among the phases r]j{c) can be seen easily from the plots of rij{c) versus c as well as 



ri{k,c) = k{k — c) versus k for a fixed value of c illustrated by Figure 5.5 Proposition 5.2 and the 
relations ( |5.4[ ) are particularly useful in order to find the asymptotic line-solitons from a given KP 
T-function as demonstrated by the example below. 

Example 5.3. Let us consider the 2x4 matrix, 

. _ /I -a 
^-[O I b c 
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Figure 5.5. The left figure shows ?7j(c) = kj{kj — c) for j = 1, ... ,4. Each 
at the intersection point of ?7i(c) — rjj{c) corresponds to c = fc^ + kj. We assume 
that the parameters kj are generic so that there is at most one intersection point 
for each c. The right figure is a plot of ?7(fc, c) as a function of k with c = fc2 + ^4, 
that is, the 77 along the dotted line in the left figure passing through the intersection 
point 772(c) = '74(c). This figure shows the order 773 < r/2 = ?74 < 771. 



where a, b, c are positive real numbers. In this case, there are six maximal minors, five of which are 
positive, namely, 

e(l,2) = l, e(l,3)=5, C(l,4)=c, e(2,4)=a, ^(3,4) = a6, 



and ^(2, 3) = 0. Then from (5.1 ) the r-function has the form, 

T =(fc2 - fci)S(l, 2) + b{k3 - ki)E{l, 3) + c(fc4 - ki)E{l, 4) 

+ a{k2 - ki)E{2, 4) + a6(fc4 - fc3)-B(3, 4) . 

Proposition |5.2| implies that the line solitons are localized along the lines x -\- cy ^constant with 
c = ki + kj = tan^'[j j]. Hence, we look for dominant exponential terms in the r-function along 
those directions. For y 3> 0, the c-values decrease as we sweep clockwise from negative to positive 
X-axis starting with the largest value c — k^ + k^. We have 771, 772 > ?73 = Vi from the order relations 
(5.4) for c = ks + k^. This means that r/i-l-r/2 is the dominant phase combination along this direction. 
Since ^(1,2) 7^ 0, T{x,y,t) (fc2 — fci)i?(l,2) implying that u w along the line [3,4], so there is 
no [3,4] line-soliton. By similar reasoning one can verify that the [1,4]- and [1, 2]-solitons are also 
impossible. Let us consider the direction c = A;2 + ^4 to check for the [2, 4]-soliton. From ( |5.4[ ) (see 
also Figure 5.5 ), 773 < 772 = 774 < 771, and since both ^(1, 2) and (^(1, 4) — a are nonzero, the r-function 
in (5.1 ) corresponds to a dominant balance of exponentials: r « (fc2 ~ ki)E{l, 2) + c(/s4 — ki)E{l, 4) 
along the line [2,4]. Therefore [2,4] corresponds to an asymptotic line-soliton as 7/ — > 00. The 
[1, 3]-soliton also exists by a similar argument. Thus, we have two asymptotic line-solitons [1,3]- 
and [2,4]-types for y ^ 0. 

We next look for the asymptotic solitons for y ^ by sweeping from the negative a:-axis to 
positive X-axis. Recall that in this case the dominant exponential E{i, j) corresponds to the least 
value of the sum rii{c) + f]j{c). It is easy to see that [1,2] and [3, 4]-solitons are impossible since 
i?(3,4) and E{1,2) are respectively, the only dominant exponentials along those directions. Then 
consider the [1, 3]-soliton. Along c = ki + k^, (5.4) implies that 772 < 771 = 773 < 7/4, and so the 
exponentials E{1, 2) and E{2, 3) would give the dominant balance. But E{2, 3) is not present in the 
above r-function because ^(2,3) — 0. So we conclude that [1, 3]-soliton does not exist as ?/ ^ 0, 
and for similar reasons, [2, 4]-soliton is also impossible. Next, checking for the [1, 4]-soliton, we have 
V2,V3 < ??i = '74 from (5.4). But as seen earlier, the dominant exponential -E(2,3) is not present in 
the r-function. However, there does exist a balance between the next dominant exponential pairs 
{£■(1, 3), i?(3, 4)} or {E{1,2), E{2,4:)} depending on whether 772 > r?3 or 772 < 773. In either case, 
there exists an asymptotic line-soliton along [1,4]. A similar argument applies along the line [2,3] 
which corresponds the other asymptotic line-soliton as y ^ 0. 
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Figure 5.6. An example of (2, 2)-soliton solution. The left left figure is at 
t = —16 and the right one at t = 16. The parameters (fci, . . . ,k4) are given by 
(—1,-0.5,0.5,2). The intermediate solitons are [3,4]-type for t = —16 and [1,2]- 
type for t — 16. Note that the middle triangular section at t = — 16 corresponds 
to the region with the dominant exponential £^(2,4), and that all five non-zero 
exponential terms in the r-function appear in the figure at t = 16. 



In summary, the r-function corresponding to the A-matrix given above, generates a KP solution 
with asymptotic hne-solitons [1,3] and [2,4] as y 0, and asymptotic line-solitons [1,4] and [2,3] 
as y ^ 0. This line-soliton solution with the parameters a = 6 = c = 1 in the A-matrix is shown in 
Figure |5.6[ 



We note that the line-solitons associated with the resonant (1, 2)- and (2, l)-soliton solutions can 
be determined in the same way as the above example by applying for the dominant balance conditions 



given by Proposition 5.2 and (5.4). We now proceed to discuss a more general characterization of all 



line-soliton solutions of the KP equation whose r- functions are given in the Wronskian form (3.101 



5.2. Characterization of the line-solitons. It should be clear from the above examples that 



a dominant exponential term determined by the relations (5.4) is actually present in the given r- 
function if its coefficient term given by a maximal minor of the A-matrix is non-zero. Thus, in order 
to obtain a complete characterization of the asymptotic line solitons, it is necessary to consider the 
structure of the N x M coefficient A-matrix in some detail. We consider the matrix A to be in 
RREF, and we will also assume that A is irreducible as defined below: 

Definition 5.4. An N x M matrix A is irreducible if each column of A contains at least one nonzero 
element, or each row contains at least one nonzero element other than the pivot once A is in RREF. 

If an iV X M matrix A is not irreducible, then the corresponding r-function gives the same KP 
solution u which is obtained from another r-function associated with a smaller size matrix A derived 
from A. One can notice from the determinant expansion in ( |5.1[ ) that 

(a) if the m-th column of A has only zero elements, then ^(mi, . . . ^tun) = if ruk = rn for 
some k, that is, the exponential E„i will never appear in the r-function; in terms of the 
chord diagram, this corresponds to a loop in the lower part of the diagram (m is a non-pivot 
index) , 

(b) if the n-th row of A has the pivot as the only non-zero element, then all ^{mi, . . . , mjv) 7^ 
contains the index n, that is, the exponential En can be factored out from the r-function; 
in terms of the chord diagram, this corresponds to a loop in the upper part of the diagram. 

So the irreducibility implies that we consider only derangements (i.e. no fixed points) of the permu- 
tation. 

We now present a classification scheme of the line-soliton solutions by identifying the asymptotic 
line-solitons as j/ — !• ±00. We denote a line-soliton solution by (7V_, 7V+)-soliton whose asymptotic 
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Figure 5.7. (iV_, iV+)-soliton solution. The asymptotic line-solitons are de- 
noted by their index pairs [e„,jn] and [inngm]- The sets {ei, 62, . . . , cat} and 
{giT 92, ■ ■ ■ , Sm-n} indicate pivot and non-pivot indices, respectively. Here — 
M — N and = N for the r-function on Gr(Af, M), and E{-, ...,•) represents the 
dominant exponential in that region. 

form consists of N_ linc-solitons as y — > — 00 and Nj^ linc-solitons for ?/ — > cxd in the xy-plane as 
shown in Figure |5.7[ The next Proposition provides a general result characterizing the asymptotic 
line-solitons of the (iV_, iV+)-soliton solutions (the proof can be found in [5]): 

Proposition 5.3. Let {ei, 62, . . . , bat} and {gi, g2, ■ ■ ■ , gM~N} denote respectively, the pivot and 
non-pivot indices associated with an irreducible, N x M , TNN A-matrix. Then the soliton solution 



obtained from the T-function in (5.1) with this A-matrix has the following structure: 

(a) For y ^ 0, there are N asymptotic line-solitons of [en, jn]-type for some jn- 

(b) For J/ <C 0, there are {M — N) asymptotic line-solitons of [im, gm]-type for some in 



An important consequence of Proposition 5.3 is that it defines the pairing map tt : [M] — > [AI] on 
i ir 

(5.5) 



the integer set [M] :— {1,2,..., M} according to 

= j„, 71 = 1,2, . . . ,iV, 
=im, m=l,2,...,M - N . 

Recall that {e„}^^j^ and {gm}m=i^ respectively, the pivot and non-pivot indices of the A-matrix 




and form a disjoint partition of [M]. Then the unique index pairings in Proposition 5.3 imply 
that the map tt is a permutation of M indices. More precisely, tt g Sm where Sm is the group 
of permutations of the index set [Af]. Furthermore, since 7r(e„) = jn > e„, n ~ 1,...,N and 



""(ffm) = *m < 5m) '71 = 1 , . . . , M — N , TT defined by ( 5.5 1 is a permutation with no fixed point 



i.e. derangements. Yet another feature of tt is that it has exactly N excedances defined as follows: 



an element I G [M] is an excedance of tt if 7r(Z) > I. The excedance set of tt in (5.5) is the set 
of pivot indices {ei, 62, . . . , ejv}- The above results can be summarized to deduce the following 
characterization for the line-soliton solution of the KP equation [8 . 

Theorem 5.5. Let A be an N x M , TNN, irreducible matrix which corresponds to a point in the 



non-negative Grassmannian Gr'^{N,M) C Gr{N,M). Then the r-function [5.1] associated with 



this A-matrix generates an [M — N, N) -soliton solutions. The M asymptotic line-solitons associated 



with each of these solutions can be identified via a pairing map tt defined by (5.5). The map tt G Sm 
is a derangement of the index set [M] with N excedances given by the pivot indices {ei, 62, ... , e^v} 
of the A-matrix in RREF. 
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Figure 5.8. An example of (3, 3)-soliton solution. The permutation of this 
solution is TT = (451263). The fc-parameters are chosen as (fci, fc2, . . . , fcg) = 
(— 1, — 0, i, 1, |). The dominant exponential for x <C is i?(l,2,5), and each 
dominant exponential is obtained through the derangement representing the so- 
lution, e.g. after crossing [5, 6]-soliton in the clockwise direction, the dominant 
exponential becomes £'(1,2,6). That is, [5, 6]-soliton is given by the balance of 
those two exponentials. 



As explained in Section |4] the derangements tt G Sm are represented by the chord diagrams with 
the arrows above the line pointing from e„ to j„ for n = 1, 2, . . . , iV, while arrows below the line 



point from gm to im for m = 1, 2, . . . , M — A^. Figure 5.8 illustrates the time evolution of an example 
of (3, 3)-soliton solution. The chord diagram shows all asymptotic line-solitons for y — >■ ±oo. 

Theorem |5.5| provides a unique parametrization of each TNN Grassmannian cell in terms of the 



derangement of Sm- This agree with the result obtained by Postnikov et al in [38ll35]. One should. 



however, note that Theorem 5.5 does not give us the indices j„ and im in the [e„,j„] and {im,9m\ 
line-solitons. The specific conditions that an index pair identifies an asymptotic line-soliton 
are obtained by identifying the dominant exponential in each domain in the xy-plane. The example 
below illustrates how to apply Theorem |5.5[ and identify all the asymptotic line-solitons for a given 
irreducible TNN ^-matrix. 

Example 5.6. Let us consider the 3x5 matrix. 




with ad — be — 0, 

where a, 6, c, d and e are positive constants, that is, the v4-matrix marks a point on Gr+(3, 5). Then 



the purpose is to find asymptotic line-solitons generated by the r-function (3.10) associated with 
this A- matrix. From Proposition |5.3[ one can see that the r-function with this matrix will produce 
a (3, 2)-soliton solution since A^ = 3 and M = 5. Moreover, the asymptotic line-solitons for this 
solution are labeled by [1, ji], [2,^2] and [4, js] for y for some ji > 1,^2 > 2 and > 4. Similarly, 
the line-solitons for y ^0 are labeled by [ii, 3] and [z2, 5] for some ii < 3 and 12 < 5. The basic idea 



to determine those indices ji,i2, is and ii, 12 is to apply Proposition 5.3 and the dominant relations 
(pl). 
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Figure 5.9. Example of shallow water waves. The upper photos are taken at a 
beach in Mexico. The lower figures show the evolution of the corresponding exact 
(3, 2)-soliton solution of (24153)-type shown in the chord diagram (see Example 
5.6 1: the left one at t = 1.5 and the right one at t = 10. (Photographs by courtesy 



of Mark J. Ablowitz.) 



Starting with the last pivot 63 = 4, it is immediate to 
We now take the next pivot 62 = 2 and find 
4, we need to check only 



Let us first consider the case for y 
find js = 5, because of js > 4 (just Proposition 5.3) 
the index j2 ■ Since the index 5 is already taken as the pair index of 63 
the cases [2,4] and [2,3]. For the existence of [2, 4]-soliton, the dominant relation (5.4) requires that 



both ^(1,2,5) and ^(1,4,5) are not zero. Calculating those minors for our A-matrix, we have 

C(l,2,5) = e^0, f(l,4,5)=d^0, 

and hence [2, 4]-soliton exists. Now we consider the case with ei = 1, that is, we have only [1,2] 
and [1,3] possibility. In the case of [1,3], we use again the dominant relation (5.4), and check 



the minors 1^(1,4,5) and ^(3,4,5) which corresponds to the dominant exponentials. We then find 
^(1,4,5) — d ^ but ^(3,4,5) = be — ad = Q. This implies that [1, 3]-soliton is impossible for 
y ^ 0. So the last one is [1, 2]-type, which can be confirmed by the condition ^(1, 4, 5) = d 7^ and 
C(2,4,5) = &^0. 

Now we consider the case for ?/ <C 0. Theorem |5.5| tells us that for the non-pivot index gi = 3, 
only the pair [1, 3] is possible (the index 2 is already taken because [l,2]-soliton exists, i.e. 7r(l) = 2). 
Then the final soliton must be [3, 5]-type from the non-pivot index (72 = 5. The last one can be 
confirmed by the least condition in ( 5.4 ) with ^(2, 3, 4) = a 7^ and ^(2, 4, 5) = 6 7^ 0. 



Thus we have a (2, 3)-soliton solution of tt = (24153)-type for the r-function (3.10) with the 



A-matrix considered. The photos in Figure |5.9| show some interacting shallow water waves, which 
we think a realization of this example. We demonstrate an exact solution whose parameters are 
given by {ki, ^2, . . . , fcs) = (—2, —1, 0, 0.5, 2) and the A-matrix with (a, b, c, d, e) — (1, 2, 1, 2, 1). 

6. (2, 2)-SOLITON SOLUTIONS 

Here we give a summary of all soliton solutions of the KP equation generated by the 2x4 



irreducible, TNN A-matrices. Proposition 5.3 implies that each of the soliton solutions consists of 
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Figure 6.1. The chord diagrams for seven different types of (2, 2)-sohton solutions. 
Each diagram corresponds to a totally non-negative Grassmannian cell in Gr(2,4). 



two asymptotic line-solitons as y ±oo. That is, they are (2, 2)-soliton solutions. We outline below 
the classification scheme for the (2, 2)-soliton solutions, and discuss some of the exact solutions in 
details for the applications discussed in the following sections. First note that for 2 x 4 matrices, 
there are only two types given by 

'l -c -d\ fl a -c^ 

.0 1 a 5 ) ^'""^ lo 1 & 



The fact that A is TNN implies that the constants a, b, c and d must be non-negative. For the 
first type, one can easily see that od = is impossible because then either ^(3,4) < or A is not 
irreducible. Then there are 5 possible cases with ad ^ 0, namely, 

(1) ad-boO, (2) ad-bc = 0, (3) 6 = , c 7^ , 

(4) c = 0,6 7^0, (5) 6 = c = 0. 

For the second type, ab due to irreducibility. Hence, we have only two cases: 

(6) CT^O, (7) c = 0. 

Thus we have total seven different types of v4-matrices, and using Theorem |5.5( we can show that 
each A-matrix gives a different (2, 2)-soliton solution which can be enumerated according to the 
seven derangements of the index set [4] — {1,2,3,4} with two cxccdances. Namely, for those cases 
from (1) to (7) we have 

(1) TT = (3412) , (2) TT = (2413) , (3) tt = (4312) , (4) tt (3421) , 
(5) TT = (4321) , (6) TT = (3142) , (7) tt = (2143) . 

In Figure |6.1[ we show the chord diagrams for all those seven cases. One should note that any 
derangement of 54 with exactly two excedances should be one of the graphs. This uniqueness in the 
general case has been used to count the number of totally non-negative Grassmann cells [351 HS] • 
Let us now summarize the results for all those seven cases of the (2, 2)-soliton solutions: 
(1) TT — (3412): This case corresponds to the T-type 2-soliton solution which was first obtained 
as the solution of the Toda lattice hierarchy This is why we call it "T-type" (see also 
[22]). The asymptotic line-solitons are [1,3]- and [2,4]-types for \y\ — >■ 00. The A-matrix is 
given by 

'1 -c -d^ 
^0 1 a b 

where a,b,c,d > are free parameters with ad — be > 0. This is the generic solution on 
the maximum dimensional cell of Gr+(2, 4), and the corresponding line-soliton has the most 
complicated pattern due to the fully resonant interactions among all line-solitons. 



A = 
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TT = (2413): The asymptotic linc-solitons are given by [1,2]- and [2, 4]-solitons for y ^ 0, 
and [1,3]- and [3, 4]-solitons for y ^ 0. The ^-matrix is given by 



A = 



1 -c 
Ola b 



where a,b,c,d > with ^(3, 4) = ad — be = 0. Note the change of the solution structure by 
imposing just one constraint ^(3,4) = to the previous case (1). 

TT — (4312): The asymptotic hnc-soUtons for this case are [1, 4]- and [2, 3]-soUtons for y ^ 0, 
and [1,3]- and [2, 4]-sohtons for y ^ 0. The ^-matrix is given by 

fl -b c 
\Q 1 a 

where a, 6, c > are free parameters. Notice that two line-solitons for y <C are the same 



as in the T-type solution (see the crossing in the lower chords in Figure 6.1). 

TT = (3421): The asymptotic line-solitons are given by [1,3] and [2,4] for y 3> 0, and for 

y ^ 0, these are the [1, 4]- and [2, 3]-solitons. The v4-matrix is given by 



A = 



1 -c 
Q I a b 



where a,b,c > are positive free parameters. This solution can be considered as a dual of 
the previous case (3), that is, two sets of line-solitons for y ^ and y ^ are exchanged 



(also notice the duality in the chord diagrams in Figure 6.1). The example discussed after 



Proposition 5.3 corresponds to this solution (see Figure 5.6 1. 

TT = (4321): The solution in this case is called the P-type 2-soliton solution which has 
asymptotic line-solitons of [1,4]- and [2,3]-types as \y\ — oo. This type of solutions fits 
better with the physical assumption of quasi-two dimensionality with weak y-dependence 
underlying the derivation of the KP equation. This is why we call it "P-type" (see p2]). 
The ^-matrix is given by 

'10 -6^ 
,0 1 a 



A 



The chord diagram indicates that those two line-solitons must have the different amplitudes, 
i.e. ^[1,4] > yl[2,3], but they can propagate in the same direction, which correspond to the 
two soliton solution of the KdV equation. 

TT = (3142): The asymptotic line-solitons are given by [1,3]- and [3, 4]-solitons for y ^ 0, 
and [1,2]- and [2, 4]-solitons for y ^ 0. The ^-matrix is given by 



A 



1 a -c 
16 



where a,b,c > 0. This solution is dual to the case (2) in the sense that the two sets of 
asymptotic line-solitons for y ^ and y <C are switched, as well as the missing minors 
are switched by C(3,4) O ^(1,2). Also note the duality between the corresponding chord 
diagrams. 

TT — (2143): This case is called the 0-type 2-soliton solution. The asymptotic line-solitons 
are of [1, 2]- and [3, 4]-types as \y\ oo. The letter "O" for this type is due to the fact that 
this solution was originally found to describe the two-soliton solution of the KP equation 
(see for example [IJ). The A-matrix for the 0-type 2-soliton solution is given by 

/I a 
yO 1 b 

Notice that this ^4- matrix is obtained as a limit c ^ in the previous one of the case (6) , 
i.e. (3142)-soliton solution. 
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-30 -20 -10 10 20 30 -30 -20 -10 10 20 30 -30 -20 -10 10 20 30 



t = -50 t = t = 50 

Figure 6.2. The time evolution of an 0-type soliton solution. Each E{i,j) in- 
dicates the dominant exponential in that region. The parameters are chosen as 

^[1,2] = ^[3,4] = 0.1 and *[3^4] = -*[i,2] = 30°. 



Now let us describe the details of some of the (2, 2)-soliton solutions, which will be important 
for an application of those solutions to shallow water problem discussed in the next Section. In 
particular, we explain how the A-matrix uniquely determines the structure of the corresponding 
soliton solution such as the location of the solitons and their phase shifts. 

6.1. O-type soliton solutions. This is the original two-soliton solution, and the solutions corre- 
spond to the chord diagram of tt = (2143). A solution of this type consists of two full line-solitons 



of [1,2] and [3,4] (see Figure 6.2 1. Note here that they have phase shifts due to their collision. Let 



us describe explicitly the structure of the solution of this type: The r-function defined in (5.1) for 
this case is given by 

r = E{1, 3) + bE{l, 4) + aE{2, 3) + abE{2, 4) , 

where a, 6 > are the free parameters given in the ^-matrix listed above. As we will show that 
those two parameters can be used to fix the locations of those solitons, that is, they are determined 
by the asymptotic data of the solution for large \y\. 

For the later application of the solution, we assume that [1, 2]-soliton has a "negative" y-componcnt 
in the wave-vector (i.e. tan^[i_2] < 0), and [3, 4]-soliton has a "positive" y-component, (i.e. 



tan^j3 4] > 0, see Figure 3.1). Then for the region with large positive x, we have [1, 2]-soliton 
in y > and [3, 4]-soliton in y < 0. 

For [1, 2]-soliton in x > (and y ^ 0), we have the dominant balance between E{1, 4) and E{2, 4). 
Then the T-function can be written in the following form, 

T w bE{l,A) +abE{2A) 

= 2be'>'+h''+'-^ cosh i (01 ^92 + e+) , 

which leads to the [1, 2]-soliton solution in the region near di w 02 for x 0, 

M = 292 inr « i(/c2 - kif sech^ ^ {Oi -~ 62 + 6I+ ) ■ 

Here the shift 9^2 (+ indicates a; > 0) is related to the parameter a in the ^-matrix (see below). 
For [3, 4]-soliton in a; > (and j/ <C 0), from the balance r w aE{2, 3) + abE{2, 4), we have 

U « i(fc4 - fc3)'sech2i(03 - 04 + 03+4). 

The shifts 9i2 and 9^^ are related to the parameters in the A-matrix, 

kA ~ kl n+ _ krt — ki „ + 
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Figure 6.3. 0-type interaction for two equal amplitude solitons. The parameters 
fci's are taken to be (fci, k2, k^, /C4) = (—1 — 10^'*, —10^^, 10"'', 1 + 10~^), which give 
^[1,2] = ^[3,4] = 5 andtan*[3^4] = -tan^'[i^2] = 1+2x10"^ (i.e. *[3_4] sa 45.0057). 
The constants a, b in the ^-matrix are chosen so that the center of interaction 
point is located at the origin, and Wmax = m(0, 0,0) ~ 1.96 and the phase shift 

^^[1,2] = ^[3,4] ~ 7.8. 



Thus the parameters in the ^-matrix can be determined by the asymptotic data of the locations of 
those [1,2]- and [3, 4] -solitons for cc » and \y\ 0. 

The most important feature of the 0-type solution is the phase shift due to the interaction of 
those two oblique line-solitons. The phase shift for [i, j]-soliton is defined by 6*^- = 6*"- — 6*^ where ± 
indicate the values for x — >■ ±00. The values O12 and O34 turn out to be the same (see for example 

where we note 



712 



InAr 



^ _ (fcs - fc2)(fc4 - fcl) - fel)(fc4 - h) ^ ^ 

°' (fc4 ~ fc2)(fc3 - fcl) {ki - k2)ik3 - ki) 

This implies that 612 = ^34 > 0, and each [i, j]-soliton shifts in x with 
(6.2) Ax[,^,] ^ " 



The positive phase shifts Aa:[i_2] > and Ax[3^4] > indicate an attractive force in the interaction. 
Figure 6.3 illustrates an 0-type interaction of two solitons which have the same amplitude, ^[1,2] = 
A [3 4] = ^, and are symmetric with respect to the y-axis, ^'[3,4] = —^[1.2] ~ 45°. Since the solution 
is close to the resonance, we have the large phase shifts Aa;jx,2] = [3 4] ~ 7.8 and the maximum 
value of the soltion M,„ax ~ 1-96 (almost four times larger than A[4_2])- 

0-type soliton solution has a steady X-shape with phase shifts in both line-solitons. One can also 
find the formula of the maximum amplitude which occurs at the center of intersection point (center 
of the X-shape), which is given by 

1-x/A^ 




(6.3) Minax = 

(see for example [HHHllllj.) Since < Aq < 1, we have the bound 



^[1,2] +^[3,4] < ""max < (^y^[L2] + y ^[3,4] 

It is also interesting to note that the formula Aq has critical cases at the values fci — k2 or 
^2 = fcs or /c3 = k^. For the case with fci — k2 or k^ = k^ (i.e. Aq = 1), one can see that one of the 
line-soliton becomes small, and the limit consists of just one-soliton solution. On the other hand, for 
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the case ^2 — (i.e. Aq — 0), the r-function has only three terms, which corresponds to a solution 
showing a Y-shape interaction (i.e. the phase shift becomes infinity and the middle portion of the 
interaction stretches to infinity). This limit has been discussed in [5S1|32] as a resonant interaction of 
three waves to make Y-shape soliton. This limit gives a critical angle between those solitons which 
can be found as follows: First let us express each kj parameter in terms of the amplitude and the 
slope, 

ki,2 = I (tan«'[i,2] T , 



- (tan'i' 
2 V 



3,4 = ^ I tan VI/ [3,4] T \/2A[3_4] 



where the angle ^'[ij] is measured in the counterclockwise direction from the y-axis (see Figure 6.3 1 
In particular, we have 



tan^'[i,2] = -^2A[i,2] +2fc2, tan^'[3,4] = ^2^(3,4] + 2/03. 

For simplicity, let us consider the special case when both solitons are of equal amplitude and 
symmetric with respect to the y-axis i.e., ^[1,2] = ^[3.4] = ^0 and ^I'p 4] — — vE'[i,2] = > 0. This 
corresponds to setting ki = —k4 and k2 = —k^. Then, for fixed amplitude Aq, the angle vJ/q has a 
lower bound given by 

tan = \/2Aq + 2k^ > \/2A^ tan ■ 
The lower bound is achieved in the limit k2 — = 0, and the critical angle is given by 

(6.4) *c = tan"^ V2A^. 

In [28], Miles introduced the following parameter to describe the interaction properties for 0-type 
solution, 

tan^-o tan^-o 

(6.5) K :— 



tan^I'^ 



With this parameter, the maximum amplitude of (6.3) for this symmetric case is given by 



(6.6) Minax = :,-^^-7=, with Ao = l-^. 

Thus, at the critical angle ^E'o = ^'c (i.e., ^ — 1), we have ttmax = 4j4o and the phase shift 6'[42] 00, 
leading to the resonant Y-shape interaction (see also [351 IH] ) . 

One should note that if we use the form of the 0-type solution even beyond the critical angle, i.e. 
^3 < ^2, then the solution becomes singular (note that the sign of E{2, 3) changes). In earlier works, 
this was considered to be an obstacle for using the KP equation to describe an interaction of two line- 
solitons with a smaller angle. On the contrary, the KP equation should give a better approximation 
to describe oblique interactions of solitons with smaller angles. Thus one should expect to have 
explicit solutions of the KP equation describing such phenomena. It turns out that the new types 
of (2, 2)-soliton solutions discussed above can indeed serve as good models for describing line-soliton 
interactions of solitons with small angles. We will show in Section |8] how these solutions are related 
to the Mach reflections in shallow water waves. 

6.2. (3142)-type soliton solutions. We consider a solution of this type which consists of two line- 
solitons for large positive x and two other line-solitons for large negative x. We then assume that 
the slopes of two solitons in each region have opposite signs, i.e. one in y > and other in y < (see 



Figure 6.4). The line-solitons for the (3142)-type solution are determined from the balance between 



two appropriate exponential terms in its r-function which has the form, 

r ^ E{1, 3) + bE{l, 4) + aE{2, 3) + abE{2, 4) + cE{i, 4) 
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Figure 6.4. The time evolution of a (3142)-type soliton solution. The line-solitons 
have ^[1,3] = ^[2,4] = 0.5 and 4* [2,4] = — ^'[1,3] = 25°. The critical angle is given by 

[3 4] — — ^'[1.2]- The parameters 



1,3] 

tan" 



45°, which also gives ^c = ^ 



in the ^-matrix are given by (6.9 1 with — = and s = 1. The amplitude of 
the intermediate [1, 4]-soliton approaches asymptotically to 1.075. 



The solution contains three free parameters a, h and c, which can be used to determine the locations 
of three (out of four) asymptotic line-solitons (e.g. two in a; ^ and one in a; <C 0). Thus, the 
parameters are completely determined from the asymptotic data on large \y\. 

Let us first consider the line-solitons in x ^ 0: There are two line-solitons which are [1, 3]-soliton 
in y S> and [2, 4] -soliton in y ^ 0. The [1, 3] -soliton is obtained by the balance between the 
exponential terms 6i?(l,4) and ci?(3,4), and the [2, 4]-soliton is by the balance between ai?(2,3) 
and ci?(3,4). Consequently, the phase shifts of [1,3]- and [2, 4]-solitons for x ^ are given by 



(6.7) d{n = In — + In - , = In -; r- + In ■ 

k^-k-^ c '^^ ki- k^ 



Now we consider the line-solitons in x <C 0: They are [3, 4]-soliton in y ^ and [1, 2]-soliton in 
?/ <C 0. The phase shifts are given respectively by 

97,, — In — ^ — In b. OZ, ~ In — ^ — In a 

fc4-fci ' A:3-fe 

We then define the parameter s (representing the total phase shifts + Q'^^ — + 6*22 ), 
(6.8) s exp (-0+ - ^34) , 



which leads to 



fci 0+ , ks- ki g+ ks - ki 



(6.9) a=-^-^se''-, fe^^^^e^ia, c=-^-^s. 

fe3 — fe2 K4 — ki K4 — K3 

The s-parameter represents the relative locations of the intersection po int o f the [1,3]- and [3,4]- 
solitons with the a;-axis, in particular, 6^^+6^^ = when s = 1 (see Figure 6.5 1. Thus the parameters 



a, 6 and c are related to the locations of [1, 3]-soliton (with ^43), of [2, 4]-soliton (with 6'Jj), and the 
intersection point of [1,3]- and [3, 4]-solitons (with s). 

Now we consider the case where the [1,3]- and [2, 4]-solitons have the same amplitude (A[i 3] = 
^[2,4] — ^0) and they are symmetric with respect to the a:-axis (^'[2.4] = ^^[1.3] = ^o)- Then in 
terms of the fc-parameters, we have 

^3 — fci = ^4 — ^2 = V^Ai- 
Also the symmetry of the wave-vectors, i.e. 4* [2.4] = ^0 = ^^[1,3] i gives 

^2 + fc4 = -(fci + fc3) = tanvJ/Q. 
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s = 10 





s = 10- 



Figure 6.5. (3142)-type soliton solution with the s-parameter. The hne-sohtons 
are given by ^[1,3] = ^[2,4] = 0.5 and VE'[2,4] = —^[1,3] = 25° (these then gives the 



other two sohtons uniquely). The parameters in the A- matrix are chosen as (6.9 1 
0. Then at s 



with 6*+ 3j 



1+ 

[2,4] 



1, all the solitons meet at the origin, i.e. the 



s-parameter shifts [1,2]- and [3, 4]-solitons. 



This implies that we have 
(6.10) 



k4 = -ki > 0, fca = ~k2 > 0. 



The angle ^'o takes the value in (0, '^c), where the critical angle is given by the condition k2 — k^ — 0, 
i.e. 

= tan" 



Notice that this formula is the same as that of the 0-type soliton solution (see (6.4)), and the 



(3142)-type exists when the K-parameter is less than one, i.e. for (3142)-type, we have 



tanvj/g 

K = , . < 1. 



From (6.10), one can easily deduce the following facts for [1, 2]- and [3, 4] -solitons in x < 0: 



(a) Those solitons have the same amplitude, i.e. 

-4[l,2] = ^[3,4] = ^(^4 - ^3)^ = ^(^4 



fc2)^ = ^ tan^ ^'o = k^Aq. 



Thus, if the [1,3]- and [2, 4] -solitons in x > are close to the y-axis (i.e. a small v^o), then 
the amplitudes of the solitons in x < are small; whereas at the critical angle '^0 = '^c, the 
solitons [1,2] and [3,4] in a; < take the maximum amplitude ^[1,2] = ^[3,4] = ^o- 
(b) The directions of the wave-vectors for the [1,2] and [3, 4] -solitons are also symmetric, i.e. 



tan^E' 



[3.4] 



tan^'i 



1,2] 



Moreover, the symmetry (6.10) implies that tan ^E* [3 4] 

*[3,4] = *r ' ' 



ki. 

- k4 



/C2 = ^2A[2A] = V^^, so 



tan 



Thus the directions of the wave- vectors for the [1, 2] and [3, 4]-solitons in a; < depend only 
on the amplitude of the solitons in a; > but not on their directions (i.e., angle of their 
V-shape). 

Let us choose the parameters in the ^-matrix for the (3142)-soliton solution appropriately, so 
that at i = all the solitons intersect at the origin (see Figure 6.4). Then for t < 0, the resonant 



interaction between [1,3]- and [3, 4]-solitons (as well as [2,4]- and [1, 2]-solitons) generates an inter- 
mediate line-soliton (called "stem" soliton) which is [1,4] soliton. The amplitude of this soliton is 
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Figure 6.6. T-type interaction with the s-parameter. The /c-parameters are chosen 
as {ki,k2,h,k4) = (-§, 



2' 2' 2 



). The A-matrix is chosen as (6.12) and (6.131 with 



^13 = ^24 = ^ ^^'^ r = 1. The s-parameter gives the phase shift for the [1,2]- and 



[3, 4]-solitons in a: < 0. 



given by 

(6.11) A[i,4] = \{k4-kif ^ ^ (v/2Ao + tan*o)' = ^o(1 + k)^ 

Note here that at the critical angle = the amplitude takes the maximum ^[1,4] = 4^o (see 

[1S1IS7]). 

For t > 0, the resonant interaction between [1,3]- and [1, 2]-solitons (as well as [2,4]- and [3,4]- 
solitons) generates an intermediate line-soliton of [2, 3]-soliton. The amplitude of [2, 3]-soliton is 
given by 



^[2,3] - 



hks - k2f = I (V^Ao - tan*o)' = ^o(l - 



2' ' 2 



Because of the symmetry (6.10 ), both [1,4]- and [2, 3]-solitons are parallel to the y-axis, i.e. tan ^'[1,4] = 
tan v]/ [2^3] = 0. 

6.3. T-type soliton solutions. There are four parameters in the A-matrix for T-type soliton 
solution. Here we explain that those parameters give the information of the locations of those line- 
solitons, the phase shift and on-set of the opening of a box. Thus three of those four parameters are 
determined by the asymptotic data on large \y\, and we need an internal data for the other one. 

Following the arguments in the previous section, one can find the phase shifts of the line-solitons 
of [1, 3] and [2, 4]: For [1, 3] -soliton in 2; > (and y ^ 0), the phase shift is calculated as 



+ - In ~ - In — 
ki — k^ b 



'13 



where D = ad ~ he = ^(3, 4). For the same soliton in a; < (and y <Si 0), we have 

a- , ^2 - fci , 
0^^ — m In c . 

So the total phase shift Oi^ 0^^ — 9^^ depends on the A-matrix unlike the cases of O- and P-types, 
and it can take any value. 

For [2, 4]-soliton in a; > (and j/ <C 0), we have 

a+ 1 ^3 - ^2 , D 
6*+ = In- — - In—. 

fc4 — fcs C 

and for the same one in x < (and y ^ 0), we have 

^24 = In 7^ T" ~ ^■ 

k4 — fci 
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Figure 6.7. T-type in terac tion with the r-parameter. 
same as those in Figure 



6.6 



The fc-parameters are the 




^24 



The A-matrix are chosen as (6.121 with 9{^^ 
and s = 1. The r-parameter gives the on-set of the box, and it does not affect the 
locations of all four line-solitons, that is, r is an internal parameter. 



Note that the total phase shift 6*24 = 0^ — 024 is the same as that for [1, 3]-soliton, i.e. the phase 
conservation along the y-axis ^13 +6*24 — ^Ta + ^24 holds. Then as in (6.8 ) for the case of (3142)-type, 
we define the s-parameter, 

s :=exp(-6i+ -6I24), 

which represents the intersection point of [1,3]- and [2, 4]-soliton. With the s-parameter, we have 

/„ „x , ^2 — ki a+ k2 — ki a+ „ ko — fcl 

(6.12) b=j^ r-■se^'^ c = -r^ — r- se'^^\ D = s. 

K4 — ki — K2 k^ — fca 

Namely, the three parameters b, c and D ~ ad — be determine the locations and the phase shift (i.e. 
the intersection point of [1,3]- and [2, 4]-solitons). One other parameter is then related to an on-set 
of a box at the intersection point (see Figure 6.7). 

In order to characterize this parameter, let us consider the intermediate solitons of [1, 4] and [2, 3]. 
First note that for t ^ 0, [1, 4]-soliton appears as the dominant balance between E{1, 2) and E{2, 4). 
Then one can find the phase shift 9^^ (here -I- indicates t > 0), 

0+ =ln^;^-ln d. 

K4 - K2 

Similarly one can get the phase shift 9^^ for i <C as 

fca - fci , D 



^14 = In 7 TT " In 

kd - k-i 



a 



Now consider the sum of Of^, i.e. 



^14 + ^14 



(fc2~fci)(fc3-fci) _^^dD 
{k^ — k2){ki — k^) a 



Also, for the [2, 3]-soliton, one can get 



q- _ , (fc2 - fci)(fc4 - k2) _ alD 
'^^-''\k,-k,){k,-k,) ''^ d ■ 



Now we introduce a parameter r in the form, 

a ki- k2 
(6.13) -=r- 



d fca — fci ' 
so that we have 



0+ + 9^i = In ^ , 9+^ + = - In (rs) . 
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Suppose that at t = 0, [1,3]- and [2, 4]-solitons in a; > are placed so that they meet at the origin, 
that is, we choose 9^^ = = 0. Also if there is no phase shifts for those sohtons, i.e. s = 1. then 
the sums become 

Oti + = In r^-{et, + 6^,). 

This impUes that at t = (and s = 1) if r = 1, then the T-type sohton solution has an exact shape 
of "X" without any opening of a box at the intersection point on the origin. Moreover, at t = if 
r > 1, then [l,4]-soliton appears in x > and [2,3]-soliton in a; < 0; whereas if < r < 1, then 



[1, 4]-soliton appears in x < and [2, 3]-soliton in x > 0. Figure 6.7 illustrates those cases with s = 1. 
The parameter r determines the exponential term that is dominant in the region inside the box. 
When r < 1, E{2, 4) is the dominant exponential term, and when r > 1 the dominant exponential is 
E{1, 3). One should note that the parameter r cannot be determined by the asymptotic data, that 
is, r is considered as an "internal" parameter. 

7. Numerical simulation and the stability of the soliton solutions 

In this section, we present some numerical simulations of the KP equation with "V-shape" initial 
wave form related to a physical situation (see for examples |3Z1 ESI US] ) ■ The main purpose of the 
numerical simulation is to study the interaction properties of line-solitons, and we will show that 
the solutions of the initial value problems with V-shape incident waves approach asymptotically 
to some of the exact soliton solutions of the KP equation discussed in the previous section. This 
implies a stability of those exact solutions under the influence of certain deformations (notice that 
the deformation in our cases are not so small). 

The initial value problem considered here is essentially an infinite energy problem in the sense 
that each line-soliton in the initial wave is supported asymptotically in either y ^ or y ^ 0, and 
the interactions occur only in a finite domain in the xy-plane. In the numerical scheme, we consider 
the rectangular domain D — {{x,y) : |x| < L^, \y\ < Ly}, and each line-soliton is matched with a 
KdV soliton at the boundaries y — iLy. The details of the numerical scheme and the results can 
be found in g^. 

We consider the initial data given in the shape of "V" with the amplitude Aq and the oblique 
angle *o > 0, 

(7.1) u{x, y, 0) = ^0 sech^ {x - \y\ tan *o) ■ 

Note here that two semi-infinite line-solitons are propagating toward each other into the positive 
x-direction, so that they interact strongly at the corner of the V-shape. At the boundaries y — iLy 
of the numerical domain, those line-solitons are patched to the KdV one-soliton solutions given by 



m(x, ±Lj^, t) = Aq sech^ \' ^ ~'~ ^ 



^ ucij.1 ' 'i'o + ^Aq. Note here that these solitons correspond to the exact one-soliton solution 
of the KdV equation with the velocity shift due to the oblique propagation of the line-soliton, 
i.e. d^u/dy^ = t&v? '^^d'^u/ dx'^ . The numerical simulations are based on a spectral method with 
window-technique similar to the method used in JIB] (see [2D] for the details). The V-shape initial 
wave was first considered by Oikawa and Tsuji (see for example [371 US]) in order to study the 
generation of freak (or rogue) waves. They noticed generations of different types of asymptotic 
solutions depending on the initial oblique angle 5*0, and found the resonant interactions which create 
localized high amplitude waves. In this section, we present the results for the cases corresponding 
to = 2 and two different angles, VE*! and \E'2 with < < ^2- where the critical angle is 
given by = tan"-'^ •\/2Ao ~ 63.4°. Then we explain these results in terms of certain (2, 2)-soliton 
solutions discussed in the previous section, and in particular, we describe the connection with the 
Mach reflection (this will be further discussed in Section [s]). 
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Figure 7.1. Initial data with V-shape wave. Each hne of the V-shape is locaUy a 
hne-sohton solution. We set those line-solitons to meet at the origin. 



The main idea here is to consider the V-shape initial wave as the part of some (2, 2)-soliton 
solutions listed in the previous section. In order to identify those soliton solutions from the V-shape 
initial wave form, let us first denote them as [ii, ji]-soliton ior y :S> and [12, J2] -soliton for y <C 0. 
Then using the relations, kj — ki = ^/2Aq = 2 and kj + fcj = tan^'o, for [i, j]-soliton and the Miles 



parameter k = tan\l/o/v2A) of (6.5), we have 
(7.2) 



ki^ — (1 ^ '^^)i ^ji — 1 ^5 

ki^ ~ (1 '^) 1 



Notice that kj^ — —ki-^ and ki^ = — fcj^ because of the symmetry in the initial wave. Moreover, at 
the critical angle 5*0 = (i-e. k = 1), we have ki^ — kj^ — 0. We also note ki^ as the smallest 
parameter and kj^ as the largest one, so that depending on the angle we obtain the following 
ordering in the fc-parameters: 

For < *o < *c (i-e. k < 1), we have 



'32 I 



implying that the corresponding chords of the and the [?2, j2]-solitons overlap. That is, [1,3] 

chord appears on the upper side of the diagram, and [2,4] chord on the lower side. This means 
that the two solitons can be identified as part of either the (3412)-type (T-type) or the (3142)-type 



solution (see the chord diagrams in Figure 6.1 1. 
For vl/^ < vj/g < I (i.e. k > 1), we have 

^ kj^ ^ ^ ki2 ^j?2 ■ 

In this case, the corresponding chords are separated, and the two solitons form part of either (2413)- 
or (2143)-type (0-type) solution. Here [1,2]- and [3, 4]-cliords appear on the upper and lower sides 
of the chord diagram, respectively. 

Then the numerical simulations show that we have the following types of the asymptotic solutions 
depending on the values ^'q: 

(a) If the angle satisfies < (i-e. k < 1), then the solution converges asymptotically to 
(3142)-type soliton solution (not T-type) 

(b) If the angle satisfies < ^'o (i.e. k > 1), then the solution converges asymptotically to an 
0-type soliton solution (not (2413)-type). 

The convergence here is in a locally defined L^-sense with the usual norm, 

1 

WfhHD) ■■= ( [[^\f{x,y)\'dxdy^' 
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where Z? C is a compact set which covers the main structure of the interactions in the solution. 
To confirm the convergence statements, we define the (relative) error function, 

(7-3) E{t) := — UoxactllL2(_Dt)/ll"oxactlli2(_Dt)i 

with the solution u*(x,y) := u{x,y,t) and an exact solution ul^^^^{x,y), where D* is the circular 
disc given by 

Dl {(x, y)€R':ix- xoit)f + (y - yo{t)f < r^} . 

The center (a;o(t), 2/o(i)) of the circular domain 13* is chosen as the intersection point of two lines de- 
termined from the corresponding exact solution. We find the exact solution ul,^^^^{x,y) by minimiz- 
ing E(t) at certain large time t — Tq: In the minimization process, we assume that the /c-parameters 



remain the same as those given by ( 7.2 1, and vary the corresponding A-matrix to adjust the solution 
pattern (recall that the A-matrix determines the locations of the line-solitons in the solution, see 
Section [5]). After minimizing E{t), that is, finding the corresponding exact solution, we check that 
E{t) further decreases for a larger time i > Tq up to a time t — Ti > Tq, just before the effects of 
the boundary enter the disc 13* (those effects include the periodic condition in x and a mismatch 
on the boundary patching). We take the radius r in 13*. large enough so that the main interaction 
area is covered for all t < Tq, but i3* should be kept away from the boundary to avoid any influence 
coming from the boundaries. The time Ti > Tq gives an optimal time to develop a pattern close to 
the corresponding exact solution, but it is also limited to avoid any disturbance from the boundaries 
for t < Ti. Thus, our convergence implies the separation of the radiations from the soliton solution, 
just like the case of the KdV equation (see the end of Section [s]). 

We also note that the convergence here implies a completion of the partial chord diagram consist- 
ing of only two chords which corresponds to the semi-infinite solitons in the initial V-shape wave. 
Namely, the asymptotic solution of the initial value problem with V-shape initial wave is given by 
an exact solution parametrized by a unique chord diagram, and the initial (partial) chord diagram 
is completed by adding two other solitons (chords) generated by the interaction. The completion 
may not be unique, and in [23|, we proposed a concept of minimal completion in the sense that 
the completed diagram has the minimum total length of the chords and the corresponding TNN 
Grassmannian cell has the minimum dimension. However, this problem is still open, and we need 
to make a precise statement of the minimal completion of partial chord diagram given by the initial 
wave profile. 

7.1. Regular reflection: k > 1. We consider the V-shape initial wave with Aq — 2 and tan = 
^,(4'o ~ 67.3°) which gives k = 1.2. Here the critical angle is *c = tan"^(2) ss 63.4°, and we 
expect asymptotically an 0-type soliton solution. The corresponding /c-parameters are obtained 



from ( 7.2 1, i.e 



(fci,fe,fc3,fc4) = (- " ' ' 



5 ' 5' 5' 5 



Figure 7.2 illustrates the result of the numerical simulation. The top figures show the direct simu- 
lation of the KP equation. The wake behind the interaction point has a large negative amplitude, 
and it disperses and decays in the negative x-direction. This shows a separation of the radiations 
from the exact solution similar to the case of KdV soliton. The steady pattern left after shedding 
the radiations can be identified as an 0-type solution. The middle figures show the corresponding 
0-type exact solution whose A-matrix is determined by minimizing the error function E{t) at t ~ 6, 

. , 1 1.91 
'0 1 0.17 



Using (6.1 ), we obtain the shift of the initial line-solitons, 

X[i,2] = a;[3_4] = -0.020. 

(Note here that because of the symmetric profile, the shifts for initial solitons are the same.) The 
negative shifts imply the slow-down of the incidence waves due to the generation of the solitons 
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Figure 7.2. Numerical simulation of V-shape initial wave for k > 1 (regular 
reflection). The initial wave consists of [1, 2]-soliton in y > and [3, 4]-soliton in 
2/ < 0, with A[i^2] = ^[3,4] = 2 and ^-q « 67.3° (*c = 63.4°. The upper figures show 
the result of the direct simulation. Notice a large wake behind the interaction point 
which extends the initial solitons. The middle figures show the corresponding exact 
solution of 0-type. The circle in these figures show the domain £)* with r = 12. 
About t = 3, the wakes seem to be out of the domain. The bottom graph shows the 
error function E(t) which is minimized at t = 6. The solid chords in the diagram 
indicates the incident solitons, and the dotted ones show the reflected solitons (i.e. 
a completion of the chord diagram |23|). 



extending the initial solitons in the negative x-direction. The phase shifts Aa;rj ,i for the 0-type 



exact solution are calculated from (6.2), and they are 



Aa:[i_2] = Aa:;[3^4] = 0.593. 

The positivity of the phase shifts is due to the attractive force between the line-solitons, and this 
explains the slow-down of the initial solitons, i.e. the small negative shifts of a^u.ji. The bottom 



graph in Figure 7.2 shows E{t) of (7.3 ), where we take r = 12 for the domain Z?* . One can see a rapid 
convergence of the solution to the 0-type exact solution with those parameters. One should however 
remark that when is close to the critical one, i.e. fc2 ~ ^3, there exists a large phase shift in the 
soliton solution, and the convergence is very slow. Note that in the limit k2 = ^3 the amplitude of 
the intermediate soliton generated at the intersection point reaches four times larger than the initial 
solitons. This large amplitude wave generation has been considered as the Mach reflection problem 
of shallow water wave [351 133 HI 123] (see also Section [S]). The chord diagram in Figure 7.2 shows 
a completion of the (partial) chord diagram: The solid chords indicate the initial solitons forming 
V-shape, and the dotted chords corresponds to the solitons generated by the interaction (see [53] 
for further discussion). 
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Figure 7.3. Numerical simulation of V-shapc initial wave for k < 1 (Mach reflec- 
tion): The initial wave consists of [1, 3]-soliton in ?/ > and [2, 4]-soliton in ?/ < 
with Aq = 2 and ^Pq — 45°. The top figures show the result of the direct simula- 
tion, and the middle figures show the corresponding exact solution of (3142)-type. 
Notice that a large amplitude intermediate soliton is generated at the intersection 
point, and it corresponds to the [1, 4]-soliton with the amplitude ^[1.4] = 4.5. The 
circles in the middle figures are D' with r = 12, which cover well the main part of 
the interaction regions up to t = 12. The bottom graph of the error function E{t) 
which is minimized at t = 10. The solid chords in the diagram indicate the incident 
solitons, and the dotted ones show the reflected solitons. 



7.2. The Mach reflection: k < 1. We consider the initial V-shape wave with Aq — 2 and — 45° 
(i.e. K = 0.5). The angle is now less than the critical angle ~ 63.4°. The asymptotic solution 
is expected to be of (3142)-type whose /c-parameters are obtained from (7.2), i.e. 



Figure 7.3 illustrates the result of the numerical simulation. The top figures show the direct 
simulation of the KP equation. We again observe a bow-shape wake behind the interaction point. 
The wake expands and decays, and then we see the appearance of new solitons which form resonant 
interactions with the initial solitons. One should note that the solution generates a large amplitude 
intermediate soliton at the interaction point, and this soliton is identified as [1, 4]-soliton with the 
amplitude ^[1,4] = 4.5. This [1, 4]-soliton is called the Mach stem in the Mach reflection [51 [53] (see 
also Section [8|. 

The middle figures in Figure |7.3[ show the corresponding exact solution of (3142)-type whose 
A-matrix is found by minimizing E(t) at t — 10, 

'1 1.92 -1.96^ 
,0 1 0.64 



A^ 
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Figure 7.4. Numerical simulation for X-shape initial wave with Aq = 2 and 
4'o = 45° (i.e. k = 0.5 < 1): The initial wave is the sum of [1, 3]- and [2, 4]-solitons. 
The top figures show the numerical simulation, and the middle figures show the 
corresponding exact solution of (3412)-type, i.e. T-type. Notice that the circle 
showing D* with r = 22 covers well the box generated by the resonant interaction 
up to i = 7. The bottom graph shows the error function E{t) of (7.3) which is 
minimized at t — 6. The four solid chords in the diagram show the asymptotic 
solitons in the incident wave, and they form a T-type soliton solution. 



Using (6.9), we obtain the phase shifts for the initial solitons of [1,3]- and [2,4]-type in a; > 0, 
and the s-parameter, 

2;[i,3] = X[2,A] = -0.01, s = 0.980. 
Those values indicate that the solution is very close to the exact solution for all the time. The 
negative value of the shifts X[ij] is due to the generation of a large amplitude soliton [1, 4]-type (i.e. 
the initial solitons slow down), and s < 1 implies that the [1, 4]-soliton is generated after t = 0. Also 
note that the [1, 4]-soliton now resonantly interact with [1, 3]- and [2, 4]-solitons to create new solitons 
[1,2]- and [3, 4]-solitons (called the reflected waves in the Mach reflection problem [551 [HI US] ) • This 
process then seems to compensate the shifts of incident waves, even though we observe a large wake 
behind the interaction point. 



The bottom graph in Figure 7.3 shows a rapid convergence of the initial wave to (3142)-type 
soliton solution with those parameters of the ^-matrix and k values given above, and for the error 
function E{t), we minimize it at t = 10 for £>* with r = 12. 

7.3. T-type interaction with X-shape initial wave. In this example, we consider an X-shape 
initial wave given by the sum of two line-solitons. For simplicity, we consider a symmetric initial 



wave with Aq = 2 and 4'o — 45° (i.e. extend the initial wave in Figure 7.3 into the negative x- 



region). Since k — 0.5 < 1, the corresponding chord diagram shows the T-type with the /c-parameters 



KP SOLITONS IN SHALLOW WATER 



41 



(/ci, ^2, fcs, A;4) = (— |,— 5, 5, §)• Although T-type sohton appears for smaller angle ^^o, one should 
not take so small value. For an example of the symmetric case, if we take ^"0 = giving twice 
higher amplitude than one-soliton case, we obtain KdV 2-soliton solution with different amplitudes 
(as can be shown by the method of 1ST). So for the case with very small angle \I/o, we expect to 
see those KdV solitons near the intersection point. However, the solitons expected from the chord 
diagram have almost the same amplitude as the incidence solitons for the case with a small angle. 
The detailed study also shows that near the intersection point for T-type solution at the time when 
all four solitons meet at this point (i.e X-shape), the solution at the intersection point has a small 
amplitude due to the repulsive force similar to the KdV solitons. Then the initial X-shape wave 
with small angle generates a large soliton at the intersection point. This then implies that our initial 
wave given by the sum of two line-solitons creates a large dispersive perturbation at the intersection 
point, and one may need to wait a long time to see the convergence. 

The top figures in Figure [7^ illustrate the numerical simulation, which clearly shows an opening of 
a resonant box as expected by the feature of T-type. The corresponding exact solution is illustrated 
in the middle figures, where the A-matrix of the solution is obtained by minimizing the error function 



E{t) of (|7J) at i 6, 

A -0.368 -0.330 
[0 1 1.198 0.123 

In the minimization process, we take a;[i 3] — a; [2, 4] due to the symmetric profile of the solution, and 



adjust the on-set of the box (see subsection 6.3). Note that the symmetry reduces the number of 
free parameters to three. We obtain 

^[1,3] = ^[2.4] — 0.025, r = 3.63, s ~ 0.350. 

The positive shifts of those [1, 3]- and [2, 4]-solitons in the wavefront indicate also the positive shift 
of the newly generated soliton of [l,4]-type at the front. This is due to the repulsive force which 
exists in the KdV type interaction as explained above, that is, the interaction part in the initial 
wave has a larger amplitude than that of the exact solution, so that this part of the solution moves 
faster than that in the exact solution. This difference may result in a shift of the location of the 
[1, 4]-soliton. The relatively large value r > 1 indicates that the onset of the box is actually much 
earlier than t = 0, and s < 1 shows the positive phase shifts as calculated from ( |6.12 1. 



The bottom graph in Figure 7.4 shows the evolution of the error function E(t) of (7.3) which is 
minimized a,t t — 6. Note here that the circular domain I?* with r = 22 covers well the main feature 
of the interaction patterns for all the time computed for t < 7. The chord diagram in the figure 
shows four asymptotic solitons in the initial wave which form a T-type soliton solution (see [20j for 
further discussion). 

8. Shallow water waves: The Mack reflection 

In this last section, we discuss a real application of the exact soliton solutions of the KP equation 
described in the previous sections to the Mach reflection phenomena in shallow water. In |28|, J. 
Miles considered an oblique interaction of two line-solitons using 0-type solutions. He observed that 
resonance occurs at the critical angle ^'c, and when the initial oblique angle ^'o is smaller than 
the 0-type solution becomes singular (recall that at the critical angle one of the exponential term 



in the r-function vanishes, see subsection 6.1). He also noticed a similarity between this resonant 



interaction and the Mach reflection found in shock wave interaction (see for example CF:48,Wh:64). 



This is illustrated by the left flgure of Figure 8.1 where an incidence wave shown by the vertical line 
is propagating to the right, and it hits a rigid wall with the angle — 5*0 measured counterclockwise 
from the axis perpendicular to the wall (see also [E]). If the angle VPo (equivalently the inclination 
angle of the wall) is large, the reflected wave behind the incidence wave has the same angle ^fg, i.e. a 
regular reflection occurs. However, if the angle is small, then an intermediate wave called the Mach 



stem appears as illustrated in Figure 8.1 The Mach stem, the incident wave and the reflected wave 



interact resonantly, and those three waves form a resonant triplet. The right panel in Figure 8.1 
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Figure 8.1. The Mach reflection. The left panel illustrates a semi- infinite line- 
soliton (incidence wave) propagating parallel to the wall with the mirror image. 
The right panel is an equivalent system to the left one when we ignore the viscous 
effect on the wall. The incident wave then forms a V-shape wave at i = as 
discussed in Section [t] The resulting wave pattern shown here is a (3142)-soliton 
solution. 



illustrates the wave propagation which is equivalent to that in the left panel, if one ignores the effect 
of viscosity on the wall (i.e. no boundary layer). At the point O, the initial wave has V-shape with 
the angle ^'o, which forms the initial data for our simulation discussed in the previous section. Then 
as we presented, the numerical simulation describes the reflection of line-soliton with an inclined 
wall, and these results explain well the Mach reflection phenomena in terms of the exact soliton 
solutions of the KP equation. 

8.1. Previous numerical results of the Boussinesq-type equations. One of the most in- 
teresting things of the Mach reflection is that the KP theory predicts an extraordinary four-fold 
amplification of the stem wave at the critical angle ^^E]. We recall the formulae of the maximum 



amplitudes which are given by (6.6) for the 0-type solution (k > 1) and (6.11) for (3142)-type 
solution (k < 1). Let a denote the amplification factor in terms of the Miles parameter n of (6.5 1, 
i.e. 

{(1 -I- k)^ for K < 1, 
, for K > 1. 
l + \/l -K-2' 

Several laboratory and numerical experiments tried to confirm the formula ( |8.1[ ), in particular, 
the four-fold amplification at the critical value k = 1 (see for example [321 I^VI Il51 [¥T1 In 
[T5] . Funakoshi made a numerical simulation of the Mach reflection problem using the system of 
equations, 

T}t + Alp + aV ■ (7]Vi)) - -A^V = 0, 

6 

^-^AV) +r; +11^12^0, 

which is equivalent to the Boussinesq-type equation ( |2.4[ ) up to this order. He considered the initial 
wave to be the KdV soliton with higher order corrections up to 0{e). In his paper, he mainly 
presented the results for the incidence waves with the amplitude = 0.05 = cio/^o ^-^id the angles 
< ^'o < f ■ He concluded that his results agree very well with the resonantly interacting solitary 
wave solution predicted by Miles. However his results on the amplification parameter a are slightly 
shifted to the lower values of the Miles K-parameter. Tanaka in [H] then re-examined Funakoshi's 
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K 

Figure 8.2. Numerical results of the amplification factor a versus the K-parameter. 
The circles show Funakoshi's result [T3], the squares show Tanaka's result [H]. The 
black dots shows the experimental results by Yeh et al 46J. 



results for higher amplitude incidence waves with a; = 0.3 using the high-order spectral method. 
He noted that the effect of large amplitude tends to prevent the Mach reflection to occur, and all 
the parameters such as the critical angle are shifted toward the values corresponding to the 
regular reflection (i.e. 0-type). For example, he obtained the maximum amplification a = 2.897 at 
K = 0.695. 

However, we claim in our recent paper [l^ that those previous results did not properly interpret 
their comparisons with the theory, and in fact their results are in good agreement with the the 
predictions given by the KP theory except for the cases near k — 1. One should emphasize that the 
KP equation is derived under the assumptions of quasi-two dimensionality and weak nonlinearity. 
Thus the key ingredient is to include higher order corrections to those assumptions when we compare 
the numerical or experimental results with the theory. In particular, the quasi-two dimensionality 
can be corrected by comparing the KP soliton with the KdV soliton in the propagation direction as 
mentioned in Section [3] More precisely, we have the amplitude correction (3.6), i.e. 

. _ ap _ 2/toAo 

°° ~ coss^-o ^ 3cos2^'o' 
where ag is the amplitude observed in the numerical computation of the Boussinesq-type equation 
(which has rotational symmetry in M^). This then suggests that the K-parameter should be evaluated 
by the following formula using the experimental amplitude uq, 

tan ^0 tan 



•2) 



^3(ao//io) cos 



Because of the quasi-two dimensional approximation, i.e. |^'o| <C 1, Miles in his paper [28] replaced 
tan^ig by V^o, and then in [TSlllI], the authors continued on to use this replacement. Then their 
computations with rather large values of gave significant shifts of the K-parameter. We then 
re-evaluate their results with our formula (8.2), and the new results are shown in Figure 8.2 Since 
Funakoshi's simulations are based on small amplitude incidence waves, his results agree quite well 
with the KP predictions. Tanaka's results are also in good agreement with the KP theory except for 
the cases near the critical angle (i.e. k = 1), where the amplification parameter a gets close to 3. This 
region clearly violates the assumption of the weak nonlinearity. Although one needs to make higher 
order corrections to weak nonlinearity, the original plots of Tanaka's are significantly improved with 
the formula (8.2). The black dots in Figure 8.2 indicate the results of recent laboratory experiments 
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Table 1. Amplification factor a for different values of k = tan 5'o/\/2Ao: 
ax=7i.i(Exp.) are the laboratory data at a; = 71.1 {x = 4.27 m), Q;a:=7i.i(KP) are 
calculated from the corresponding KP exact solutions at t — 41.05, and ax=oo(i^P) 
are from (6.6) and (6.11). In the row of Aq — 0.413 with vj/g = 30°, the values of a 
in the brackets are obtained at a; = 50.8, because of the wave breaking immediately 
after this point; hence, the greater amplification cannot be realized [iT. 



K 


Ao 


*o 


aa;=71.l(Exp.) 


a.=7i.i(KP) 


a.=oo(KP) 


1.392 


0.086 


30° 


2.10 


2.36 


2.36 


1.242 


0.108 


30° 


2.13 


2.51 


2.51 


1.017 


0.161 


30° 


2.24 


3.38 


3.38 


0.887 


0.212 


30° 


2.33 


2.43 


3.56 


0.731 


0.312 


30° 


2.52 


2.61 


2.99 


0.722 


0.127 


20° 


1.89 


1.84 


2.96 


0.635 


0.413 


30° 


(2.48) 


(2.54) 


2.67 


0.591 


0.189 


20° 


1.95 


1.93 


2.53 


0.516 


0.249 


20° 


1.99 


2.08 


2.30 


0.425 


0.367 


20° 


2.01 


1.99 


2.03 



done by Yeh and his colleagues [46j. We will discuss their experimental results in the next section. 

Before closing this section, we remark on the length of the Mach stem (i.e. the intermediate 
soliton of [1, 4]-type). From (3142)-soliton solution, one can find the point (x* , y*) of the interaction 
of the triplet with [1,3]-, [3,4]- and [1, 4]-solitons ^8' , 

X* = ^ (tan + taii-^of t ^ ^{l + K)^t, 



V* = ^(tan*c-tan«'o)t 



(1 - K)i, 



(see Figure 8.1 1. In the physical coordinates with (i*, y*, i), we have 

/ ClQ cos^ ^1^0 



Co 1 IH j^^^ (1 + k) ]t, 



y* 



Co] 



3/in 



(1 - k) i. 



The angle $ in Figure 8.1 is then given by tan<I> — which is approximated in [211 [T5] by 



tan<i> 



Using the corrected formula tan <I> 
(see Figure 8 in [H]). 



y*/i*, one can see again a good agreement with the KP theory 



8.2. Experiments. Recently, Yeh and his colleagues [IHI performed several laboratory experiments 
on the Mach reflection phenomena using 7.3 m long and 3.6 m wide wave tank with a water depth 
of 6.0 cm. Here we briefly describe their results and show that our KP theory can predict very well 
the evolution of the waves observed in the experiments. 

The wave tank is equipped with 16 axis directional-wave maker system along the 3.6 m long 
side wall marked by a; = 0. An oblique incident solitary wave is created by driving those 16 
paddles synchronously along the sidewall, and the wave maker is designed to generate a KdV soliton 
with any heights before the breaking. The temporal and spatial variations of water-surface profiles 
are measured by the Laser Induced Fluorescent (LIF) method (a highly accurate measurement 
technique). The water dyed with fluorescein (green) fluoresces when excited by the laser sheet. 
The illuminated image of the water-proflles are recorded by a high-speed and high-resolution video 
camera. 
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Figure 8.3. Two views of the temporal variation of the water-surface profile in the 
y-direction (perpendicular to the wall) at x = 71.1 [37]: The top panels show the 
experimental result, and the bottom ones show the corresponding (3142)-type exact 
soliton solution of the KP equation. The incident wave amplitude Aq = 0.212, and 
the angle 5*0 = 30° . The amplification factors obtained from the experiment and the 
exact solutions are close, and they are aa;=7i.i(Exp.) — 2.33 and Q;a:=7i.i(KP) = 2.43 
(see Table [T]). 



8.2.1. The Mach reflection. In Table [T] their experimental results of the amplification factors a are 
compared with those obtained from the exact solutions of the KP equation (i.e. 0-type for k > 1 
and (3142)-type for k < 1). The waves were measured at 5: = 4.27 m {x = x/Hq = 71.1) which is the 
farthest measuring location in the experiments, except for the case with Aq = 0.413. In the later 
case, the a values in the brackets are measured at a; = 50.8 because of the wave-breaking (notice 
that at this point a = 2.48 implies — 1.02). We calculate the corresponding KP exact solution at 
t — 41.05 (recall here that the relations ( |3.8[ ) gives x — coi = Hqx and coi = ^t). The amplification 
factor a is still growing along the propagation direction, and the values obtained from the exact 
solutions are in good agreement with the measurements. We note here that near the critical case 
(i.e. K = 1) the growth of the stem amplitude is very slow and a.t x = 71.1 (x — 4.27 m) the 
amplification factor is only achieved about 65%. Also for the cases with small oblique angle, i.e. 
^0 = 20° in Table [ij the amplification factor a grows slower when the incidence wave amplitude 
is smaller, that is, at x = 71.1, a is almost constant (slightly decreases) as k increases. However 
the asymptotic value of a for large x increases as k increases. This means that the observed waves 
are still in the transient stage, and a longer tank is necessary to observe further growth of the 
amplification factor (see [17] for the details). In Figure 8.3 we show the image of the wave-profile 
at a; = 71.1 for the case when the incident wave amplitude have Aq = 0.212 and = 30°. The 
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corresponding exact solution with those parameters is of (3142)-type (i.e. k = 0.887 < 1). The 
upper panels show two views in different angles of the temporal variation of the wave profile of the 
experiment at x = 71.1, which is made from 250 slices of the spatial profiles (100 slices per second) 
with approximately 3000 pixel resolution n the y-direction. As expected form the (3142)-type exact 
solution, the stem-wave formation is realized, in which the incident and refiected waves separate 
away from the wall by the stem-wave. The lower panels show the corresponding (3142)-type exact 
solution at t — 41.05 (the x-coordinate is converted to the ^-coordinate, using ( |3.8| )). Here the 
fc-parameters are (fci, fcz, /ca, ^4) = (-0.614,-0.037,0.037,0.614) from Aq = 0.212 and ^-q = 30°. 
Then we calculate the A-matrix using (6.9 1 with s = 1, and take 



A 



1 8.797 -1.128 
1 0.530 



This choice of the ^-matrix places the incidence wave crossing at the origin at t = 0, i.e. 6'^ gj — 
(6.9). We see a good agreement between the experiment and the KP theory. At 



12,4] 



= in 



X = 71.1, the wave-profile observed in the experiment is close to the corresponding exact solution 
of (3142)-type, that is, the radiations generated at the beginning stage dispersed and well separated 
from the main part of the wave-profile as predicted in the numerical simulation (see subsection 7.2 1. 



8.2.2. T-type interaction. Figure [8^ shows a preliminary result for T-type interaction pattern gen- 
erated in the same tank by Yeh and his collaborators. The T-type solution is the most complex and 
interesting soliton solution associated with the r-function on Gr+(2,4). The initial wave has the 
V-shape (half of the X-shape) , then other half of the X-shape is generated by the line-soliton with op- 
posite angle. The upper panels of Figure 8.4 show the evolution of the wave pattern with Aq — 0.431 



and ^'o = 25°. Behind the crossing wave form (the right side in the figure), the large wakes are 
generated at the early stage of the evolution, but they eventually separate from the main pattern 
of the T-type interaction, as we observed in the numerical simulation. The figures clearly show the 
formation of a box pattern as expected by the KP theory. The lower panels show the corresponding 
T-type exact soliton solution whose parameters are (fci, fc2, ^3, ^4) = (-0.697,-0.231,0.231,0.697) 
and the v4-matrix given by (6.12) and (6.131 with 6'^ = ^'[9 41 = 0, s = 40 and r = 1, i.e. 



A = 



1 -40.34 
1 24.07 



[2,4] 

-24.07 
13.37 



The large value of the s-parameter indicates large phase shift of the incident line-solitons (see Figure 
6.6 ), that is, the parts of line-solitons in the right side (i.e. behind the interaction point) were created 
with some delayed time. In finding those parameter values, we did not make a precise minimization 
of certain error function, like the one given in (7.3). In a future communication, we hope we will be 



able to develop a method to determine the exact solution from the experimental data, that is, the 
inverse problem of the KP equation. 
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Figure 8.4. T-type solutions generated in the water tank [4^. The experimental 
results are shown in the upper panels. Those figures are made by combining the real 
image of the wave profiles with their mirror images of the wall at y = (the center 
horizontal line). The lower panels show the corresponding T-type exact solution 
of the KP equation. The initial wave has = 2.1 cm {Aq — 0.43) and '^q = 25° 
(i.e. K = 0.503). The exact solution is plotted in the xy-plane, and the solitons are 
propagating to the left. (The experimental figures by courtesy of Harry Yeh.) 
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